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Abstract. We give the first rigorous derivation of the celebrated Ginzburg- 
Landau (GL) theory, starting from the microscopic Bardeen-Cooper-Schrieffer 
(BCS) model. Close to the critical temperature, GL arises as an effective 
theory on the macroscopic scale. The relevant scaling limit is semiclassical in 
nature, and semiclassical analysis, with minimal regularity assumptions, plays 
an important part in our proof. 



1. Introduction and Main Results 

1.1. Introduction. In 1950 Ginzburg and Landau [7j introduced a model of su- 
perconductivity that has been extremely successful and is widely used in physics, 
even beyond the theory of superconductivity. It has a rich mathematical struc- 
ture that has been studied in great detail, and has inspired the development of 
many interesting new concepts. Ginzburg and Landau arrived at their model in a 
phenomenological way, describing the macroscopic properties of a superconductor, 
without the need to understand the microscopic mechanism. 

In 1957 Bardeen, Cooper and Schrieffer ;.2; formulated the first microscopic ex- 
planation of superconductivity starting from a many-body Hamiltonian. In a ma- 
jor breakthrough they realized that this phenomenon can be described by a pairing 
mechanism. Below a critical temperature, a superconducting paired state forms due 
to an instability of the normal state in the presence of an attraction between the 
particles. In the case of a metal, an effective attraction between the electrons arises 
through an interaction with phonons, the quantized vibrations of the underlying 
lattice formed by atoms. 

A connection between the two approaches, the phenomenological GL theory and 
the microscopic BCS theory, was established in 1959 when Gorkov 8 explained 
how, close to the critical temperature, the GL theory arises from the BCS model. A 
simplified version of his explanation was later given by de Gennes in his textbook [6] . 
The purpose of this paper is to give the first fully rigorous mathematical derivation 
of GL theory from the BCS model. In our approach it is not necessary to restrict 
attention to the particular interaction considered in the original BCS model, indeed 
we allow for a large class of attractions among the particles. In particular, we 
allow for local two-body interactions, in contrast to the simpler (non-local rank 
one) interaction considered by BCS. Such local interactions are relevant for the 
description of cold gases of fermionic atoms, which are of great current interest. 
These atoms are electrically neutral and the corresponding pairing mechanism is 
relevant for superfluidity rather than superconductivity. 
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The BCS model is a considerable simplification of the full many-body problem. 
In the latter a state is described by a complicated wave function of a macroscopic 
number of variables. In the BCS model, all the information about the system is 
contained in quantities depending on only two variables; the reduced one-particle 
density matrix 7, i.e., a positive trace class operator on the one-particle space, and 
the Cooper pair (two particle) wave function a, which is non-zero only below the 
critical temperature. The GL model is yet much simpler, as it describes the system 
by a single function ip of one variable only, which satisfies a non-linear second order 
PDE, the GL equation. This function ip only describes macroscopic variations in the 
system, whereas the BCS states 7 and a exhibit both microscopic and macroscopic 
details. Hence GL theory represents a significant simplification as to the more 
complicated BCS theory. The relation of these three theories is analogous to atomic 
physics where quantum mechanics, Hartrce-Fock theory and Thomas-Fermi theory 
are models of the same type of increasing simplicity. In contrast to the atomic 
case, where the Hartree-Fock approximation is mathematically well understood, it 
remains an open problem to rigorously establish the relation of the BCS model to 
the full many-body quantum-mechanical description. 

The BCS theory may be thought of as a variational theory, which utilizes a 
special class of trial states known as quasi-free states. This point of view was em- 
phasized by Leggett [12] . The mathematical aspects of this theory in the translation 
invariant case were studied in [5j [9l [TTJ [121 [IF • In this paper we study the non- 
translation invariant case, where weak external fields are present that vary only on 
the macroscopic scale. We show that close to the critical temperature Tc, where a 
is small, the macroscopic variations of a are correctly described by the GL theory. 
More precisely, we assume that the microscopic scale of the system is of order /i <C 1 
relative to the macroscopic scale. Variations of the system on the macroscopic scale 
cause a relative change of the energy of order h^, hence the external fields should 
change the energy to the same order. This change is correctly described by the GL 
theory to leading order in h if (Tc — T)/Tc is of the order , with T the actual 
temperature. We also prove that the GL wave function correctly describes the 
macroscopic behavior of the BCS state. The BCS Cooper pair wave function in 
this parameter regime equals 

"(a;, y) = \ {-ipihx) + ilj{hy)) ao{x - y) (1.1) 

to leading order in h, where ao is the unperturbed translation invariant pair func- 
tion. Equivalently, ^{i{j{hx) + 'tp{hy)) could be replaced by 'tp{h{x + y)/2) to leading 
order. In particular, ip describes the center-of-mass motion of the BCS state. The 
role of /i as a semiclassical parameter can be understood as follows. If one rescales 
the arguments to macroscopic variables and thinks of a{x/h,y/h) as the integral 
kernel of an operator, it is the quantization of the semiclassical symbol ij}{x)aa{p) . 
It would, in fact, be the Weyl quantization if we had used the alternative center- 
of-mass representation. For technical reasons we found it more convenient to work 
with the representation (|l.ll) . 

Our approach is motivated by the work of de Gennes [5]. While de Gennes 
studied the emergence of the GL equation from the BCS equations, it is important 
for our rigorous analysis to utilize that these equations are Euler-Lagrange equa- 
tions for variational problems. We give precise bounds on the lowest energy of the 
BCS functional and the connection of the corresponding minimizer with the GL 
minimizer. The exact statement of our result is given in Subsection 11.41 
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1.2. The BCS Functional. We consider a macroscopic sample of a fermionic 
system, in d spatial dimensions, where 1 < d < 3. Let /i G M denote the chemical 
potential and T > the temperature of the sample. The fermions interact through a 
local two-body potential V. In addition, they are subject to external electric and/or 
magnetic fields. Neutral atoms would not couple to these fields, of course, but 
there can be other forces, e.g., arising from rotation, with a similar mathematical 
description. In BCS theory the state of the system can be conveniently described 
in terms of a 2 x 2 operator valued matrix 

^ (1.2) 

satisfying < T < 1 as an operator on L'^{R'^) L'^{R'^) = L'^{R'^) (g) C^. The bar 
denotes complex conjugation, i.e., a has the integral kernel a{x, y). In particular, T 
is hermitian, implying that 7 is hermitian and a is symmetric, i.e., 7(0;, y) — "f{y, x) 
and a{x,y) — a{y,x). There are no spin variables in T. The full, spin dependent 
Cooper pair wave function is the product of a with an antisymmetric spin singlet. 
This is why a itself is symmetric so that we obtain the antisymmetric fermionic 
character of the full, spin-dependent, pair wave function. 

The general form of the BCS functional for the free energy of such a system is, 
in suitable units, 

(^{-iV + A{x)f - fi + W{x)'^-f -TS{T)+ 1 V{x^y)\a{x,y)\^dxdy. 



Tr 

Here, A is the magnetic vector potential, and W is the external electric potential. 
The entropy equals 

S{r) = -Tr [rinT] , (1.3) 

where the trace is now both over and L^(R'^). The BCS state of the system is 
a minimizer of this functional over all admissible states F. 

As explained in detail in [9l Appendix A], the BCS functional can be heuristically 
derived from the full many-body Hamiltonian via two steps of simplification. First, 
one considers only quasi-free states, and second one neglects the resulting direct 
and exchange term in the interaction energy. The latter terms are considered 
unimportant in the physically relevant parameter regimes (161. 

We are interested in the effect of weak and slowly varying external fields, as al- 
ready explained in the previous subsection. More precisely, A{x) should be replaced 
by hA{hx) and W{x) by h^W{hx). In order to avoid having to introduce boundary 
conditions, we assume that the system is infinite and periodic with period in 
all d directions. In particular, A and W should be periodic. We also assume that 
the state F is periodic. The aim then is to calculate the free energy per unit volume. 

We find it convenient to do a rescaling and use macroscopic variables instead of 
the microscopic ones. The rescaled BCS functional has the form 



_^BCS 



(F) := Tr ((-i W + hA{x)f - m + h^W{x)^ -f - T S{T) 

V{h-\x~y))\a{x,y)fdxdy (1.4) 



where C denotes the unit cube [0, 1]'', and Tr stands for the trace per unit volume. 
More precisely, if _B is a periodic operator (i.e., it commutes with translations by 1 
in the d coordinate directions), then Tr J5 equals, by definition the (usual) trace of 



4 R.L. FRANK, C. HAINZL, R. SEIRINGER, AND J. P. SOLOVEJ 

X-B, with X the characteristic function of C. The location of the cube is obviously 
of no importance. It is not difficult to see that the trace per unit volume has 
the usual properties like cyclicity, and standard inequalities like Holder's inequality 
hold. This is discussed in more detail in Section [3] 

We make the following assumptions on the functions A, W and V appearing in 
(ILI. 

Assumption 1. We assume both W and A to he periodic with period 1. We further 
assume that W{p) and \A(p)\{l + \p\) are summable, with W{p) and A{p) denoting 
the Fourier coefficients of W and A, respectively. In particular, W is bounded and 
continuous and A is in C^{M.'^). 

The interaction potential V is assumed to be real-valued and reflection-symmetric, 
i.e., V{x) = V{~x), with V £ LP(R''), where p = 1 for d ^ 1, p > 1 for d ^ 2 and 
p = 3/2 for d = 3. 

Our results presumably hold under slightly weaker regularity assumptions on 
W and A. For the sake of transparency we shall not aim for the weakest possible 
conditions, but rather try to keep the proofs to a reasonable length. 

The assumption on V is the natural one to guarantee relative form-bounded- 
ness with respect to the Laplacian. A slower decay at infinity could possibly be 
accommodated. Our method also works for non-local potentials, which appear 
naturally in the theory of superconductivity. For simplicity, we are not aiming for 
the most general setting and work with Assumption [T] from now on. 

We note that a periodic magnetic field B, satisfying Maxwell's equation W-B ~ 0, 
can be described via a periodic vector potential A as i? = V A A if certain flux 
conditions are satisfied. For d = 2, B has to have zero average, while for d = 3 the 
flux through the boundaries of the unit cell has to vanish. This follows from the 
fact that the relevant de Rham cohomology (closed two-forms modulo exact ones) 
equals for the n-dimensional torus. 

1.2.1. The Translation- Invariant Case. In the translation invariant case, where 
W ^ A — 0, it makes sense to restrict J^^*-^^ to translation invariant states F. 
This is the case studied in detail in In particular, it was shown in [5] that 

there is a critical temperature Tc > such that for T > T^, the BCS functional 
jtBCS jg minimized for a = and 7 = (1 -I- exp((— /i^V^ — ijl)/T))~^, i.e., 7 is 
the one-particle density matrix of a free Fermi gas. The critical temperature Tc is 
determined by the unique value of T such that the operator 

Kt + V{x) (1.5) 

on Ll^^iW^), the reflection-symmetric square-integrable functions, has zero as its 
lowest eigenvalue0 Here, Kt denotes the operator 




^The results in [9] are worked out in three dimensions, but analogous results are easily seen to 
hold in one and two dimensions. 

^In 9 the symmetry constraint on a was not explicitly enforced. The results hold equally if 
one works only in the subspace of reflection symmetric functions in L^(R'*), however. 
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Note that Kt is monotone increasing in T and Kt > 2T. (If Ko + V = | — — ^| + 
> 0, then Tc = 0.) In particular, the essential spectrum ol Kt + V is [2T, oo), 

and hence an eigenvalue at is necessarily isolated. 

In the following, we shall assume that Tc > and that the ground state of (|1.5p 

a,t T — Tc is non-degenerate. We emphasize that Tc is independent of h. 

Assumption 2. The potential V is such that Tc > and that Kt^ + V has a 
non- degenerate ground state eigenvalue 0. 

According to |TT] this assumption is satisfied if < (and not identically zero), 
for instance. In the case that V is invariant under rotations, the non-degeneracy 
assumption means that the minimizing function has angular momentum zero. We 
denote the eigenfunction of Kt^ + V corresponding to eigenvalue zero by ao. 

1.3. The GL Functional. Let ip e i/i„(R'')= the periodic functions in H^^^iR''). 
For Bi a positive dx d matrix, B2 eM. and B3 > 0, the GL functional is defined as 



/ ((_iV + 2A(x))V) •Bi((-iV + 2.4(x))V') 



+B2Wix)\iPix)\^ + B3\l-\ijix)\^\ dx. (1.7) 



Note the coefficient 2 in front of the vector potential A. It is due to the fact that ip 
describes pairs of particles, and the charge of a pair is twice the particle charge. The 
relevant coefficients Bi, B2 and B3 will be calculated below from the BCS theory, 
see Eqs. (|1.19|) " (ll.2ip . In the rotation invariant case (i.e., for radial potentials V), 
the matrix Bi is a multiple of the identity matrix. 

The ground state energy of the GL functional will be denoted by 

= inf : V e H^eM')} ■ (1-8) 

It is not difficult to show that under our assumptions on A and W, there exists a 
corresponding minimizer, which satisfies a second order differential equation known 
as the GL equation. The mathematical aspects of the GL functional have been 
studied extensively in the literature; see [H [21] and references therein. 

1.4. Main Results. Recah the definition of the BCS functional J"^^^ in (fL4)l . 
Recall also that admissible states F are of the form (|1.2|) . are periodic with period 
1 and with < 7 < 1. We define the energy F^'^^{T,ii) as the difference between 
the infimum of T^^^ over all admissible F and the free energy of the normal state 

where 70 = (1 + e((-''''^+'''4(a;))'-(-'i'iv(3;)-M)/T)-i^ i 

F^C"(r,M) = inf -^^^'(r) - J-^^S(Fo) . (1.10) 

Note that Fq is the minimizer of the BCS functional in the absence of interactions, 
i.e., when V = 0. We have 

J-BCS(p^) = -TTr In (l + cxp (- ((-i/iV + hA{x)f - /i + h'^W(x)'^ /t)) , 

(1-11) 

which is 0{h '*) for small h. The system is said to be in a superconducting (or 
superfluid, depending on the physical interpretation) state if F^'~"^ < 0. 
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THEOREM 1. As above, let Tc > denote the critical temperature in the transla- 
tion invariant case, and assume Assumptions^^ and\^ Let D > 0. For appropriate 
coefficients Bi, B2 and B3 (given in U.19\) - !il.21\) below), we have, as h 0, 

F^'^^iTcil - h^D),fi) ^ h'^-'^ {E^^ - B3 + e) , (1.12) 

with e satisfying the bounds 

( h^/^ ford=l 

const /i> e> -CL :=-constx i h^/^ [\n{l/h)]^^^ ford = 2 (1.13) 

[ h^/^ for d = 3. 

Moreover, ifT is an approximate minimizer of T^'^^ at T — Tc(l — h^D), in the 
sense that 

_^BCS(p) < J-BCS(r^) + (£;GL _ + e) (1.14) 

for some small e > 0, then the corresponding a can be decomposed as 

a = ■^(V'(x)ao(-i/iV) + ao(-iW)-0(a;)) + cr (1.15) 

with £'^^{ip) < E'^^ + e + cl, ao the (appropriately normalized; see U.18]) ) zero- 
energy ground state of U.5\) at T = Tc, and 

\a{x,y)\'^ dxdy < const — —. (1-16) 

To appreciate the bound (|1.16l) . note that the square of the L^{C x M'') norm of 
the first term on the right side of (jl.lSp is of the order h^~'^, and hence is much 
larger than the one of a for smaU h. Theorem [1] thus justifies the claim made in 
the Introduction in (|l.ip . 

For smooth enough A and W , one could also expand ^-"^^^(ro) in (II. lip to order 
and thus obtain the total energy infr J^^'~"^(r) to this order. Our approach 
bounds directly the difference (jl.lOp . however, and hence it is not necessary to 
compute J^B'~"^(ro) in detail, and to make the corresponding additional regularity 
assumptions on A and W . 

Our proof shows that the coefficients Bi, B2 and Bj, in the GL functional are 
given as follows. Recall that ao denotes the unique ground state of Kt^ + V . It 
satisfies ao(a;) = a^i^—x), and we can take it to be real. Let t denote the Fourier 
transform of 2KT^aQ = —2Vao, i.e., 

t(p) = -2(2^)-'^/2 /■ v{x)aQ{x)e"P-'=dx. (1.17) 



We normalize it, i.e., choose the normalization of ao, such that 

where Pc = ^/Tc, D = (1 — T/Tc)/h'^ as above, and gi denotes the function 



o2z 



2ze^ - 1 



z2(l + e^)2 
Note that gi{z)/z > 0. Define also 

2e^ (e^ - 1) 

92{z) = — -3- 

z(e^ + 1) 
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(compare with p.6p ~ ()2.9p ). Then Bi is the matrix with components 

(Bi)., = § {5^MPc{q^ - m)) + 2/3cg.9, 92{Pc{q^ - . (1.19) 

Note that in case V is radial, also ag and t are radial and thus Bi is proportional 
to the identity matrix. Moreover, B2 and -B3 are given by 

B2^^f t{qr9imq'-^^))j^ (1.20) 

and 

respectively. Alternatively, could be written using the normalization (jl.lSp . In 
particular, B^/\B2\ and i?3/|jBi|| are proportional to Z?, and hence proportional to 
the difference of the temperature to the critical one. 

Note that i?3 > since gi{z)/z > for all z g M. The coefficient B2 can, in 
principle, have either sign if /i > 0, however. It has the same sign as the derivative 
of Tc with respect to /x. This is not surprising. The external potential plays the 
role of a local variation in the chemical potential. If increasing /i increases T^, the 
pairing mechanism should be enhanced at negative values of W and the density of 
Cooper pairs should be largest there. If increasing /i decreases Tc, however, the 
situation is reversed. 

To see that the matrix Bi is positive, we calculate 

2Re(ao|xi(i^T, + V)xj\ao) = (ao|[x.i, [/sTt,, a;j]]|ao) 



.9i(/3,(g^~/.)) 
sinh(/3c(p^ - Af) 



8(27r)'*Bl,,, - 2Pt / t{qY q,qj '7 ' , o \ dq ■ (1.22 



The last term on the right side defines a negative matrix, while the left side defines 
a positive matrix, since Kt^ + V >Q. Hence Bi is positive. 

1.5. Outline of the paper. The structure of the remainder of this paper is as 
follows. In Section [2] we shall state our main semiclassical estimates. These are 
a crucial input to obtain the bounds in Theorem [1] The leading terms in our 
semiclassical expansion can, in principle, be obtained from well-known formulas in 
semiclassical analysis, but the standard techniques do not apply directly because we 
are forced to work with rather minimal regularity assumptions. We shall formulate 
the main results in separate theorems; see Theorems [2] and |3] below. Their proofs, 
which are rather lengthy, will be given in Sections [8] and [9l respectively. Some 
technicalities are deferred to the appendix. 

Section [3] explains various inequalities for the trace per unit volume which will be 
used throughout the proofs. Properties of the ground state of (|1.5p at T = Tc, 
are derived in Section |4l It is shown that i, defined in (|1.17p above, is smooth and 
has a suitable decay at infinity. The results imply, in particular, that the coefficients 
prTO)) ~ pr^ are well-defined and finite. 

An upper bound on i^^^^s ^jjj derived in Section [SJ using the variational 
principle. An important input will be the semiclassical estimates of Theorems [2] 
and[3l 

Sections [6] and [7] contain the lower bound. In the first part, the structure of 
an approximate minimizer is investigated, which leads to a definition of the order 
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parameter ip. This structure is then a crucial input to the second part, where also 
the semiclassical estimates of Section [2] enter. 

Throughout the proofs, C will denote various different constants. We will some- 
times be sloppy and use C also for expressions that depend only on some fixed, 
/i-independent, quantities like Tc or ||M^||oo, for instance. 



2. Semiclassical Estimates 

One of the key ingredients in the proof of Theorem [T] is semiclassical analysis. 
Choose a periodic function tp £ iJj^^(R'^) and a sufficiently nice function t and let 
A be the operator 

A = {tp{x)t{^ih\/) + t{-ih\/)tP{x)) . (2.1) 
It has the integral kernel 

A(x, y) = -^1^ + 4'{y))t{h-\y - x)) . (2.2) 

Our convention for the Fourier transform is that f{p) — (27r)^'^/^ J^^ f{x)e~^P'^dx. 
We shall assume that 

£ 2.2p/(p-i)(Rd) (2.3) 
(with p defined in Assumption [1} and that 

^.,<oo (2.4, 

for all 7 e {0,1 ...,4}''. For simplicity, we also assume that t is reflection- 
symmetric and real- valued. For the function t in (I1.17p . these assumptions are 
satisfied, as will be shown in Section 01 
Let Ha be the operator 



Ha 



{-ihV + hA{x))- ii + h^W{x) A 

A -{ihV + hA{x)f + ^i-h^W{x) 

(2.5) 

on L^(M'*) C^, with A and W satisfying Assumption [1] In the following, we 
will investigate the trace per unit volume of functions of . Specifically, we are 
interested in the effect of the off-diagonal term A in Ha, in the semiclassical regime 
of small h. 

THEOREM 2. Let 

/(z) = -ln(l + e-^) , (2.6) 

and define 

z z 

gAz) - -9,{z) + - ,2(i + e^)2 (2-8) 

,/ N 2 , ^ /"'(z)-/"'(-z) 2e^(e^-l) 

g^(z) ^ g[{z) + - g,{z) ^ ' ^ ' ^ ^ ' = -3^. 2.9 

z z z^e'' + 1) 



an. 



d 
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Then, for any /3 > 0, the diagonal entries of the 2x2 matrix-valued operator 
/(/3iJA)^/(/5i?o) are locally trace class, and the sum of their traces per unit volume 
equals 

J Tr [/(/3i/A) - fWHo)] = h'E^ + h^E^ + 0{h^) {\mW) + UWW)) 

+ 0{h^){\\4,tHHc) + \\nH^c)) : (2-10) 

where 

El = ~^\ml I t{qf goWiq' - /i)) (2.11) 

2 Jjjd (277)" 



^2 = -f E / t{q)[d,dkt]{q)goiP{q'-f^))j^ 

+ ^ E ((^J- + '^^^Mi.dk + 2zAfc)V') 



X / t{qf {Sjkgmq^ - M)) + 2/3q,gfcg2(/3(g' - m))) 



dq 
(2^ 



i J^d (Ztt) 



d 



+ ^Mtf t(,)^ ^^W-/-)) ^. (2.12) 

The error terms in i2.10\) of order and h^ are bounded uniformly in (3 for /3 in 
compact intervals in (0,00). They depend on t only via hounds on the expressions 
KM and |0). 



Here, we use the short-hand notation HV'Ilp for the norm on LP{C). Likewise, 
( • I • ) denotes the inner product on L^{C). 

In general, the operator f{f3H/s) — f{(3HQ) wiU not be trace class under our 
assumptions on t and ip. Hence the trace in (|2.10l) has to be suitably understood 
as the sum of the traces of the diagonal entries. This issue is further discussed in 
Section [5j 

The expressions Ei and E2 are the first two non- vanishing terms in a semiclassical 
expansion of the left side of (|2.10l) . They can be obtained, in principle, from well- 
known formulas in semiclassical analysis [UJ [5D] . The standard techniques are not 
directly applicable in our case, however. This has to do, on the one hand, with our 
rather minimal regularity assumptions on W, A, tp and t and, on the other hand, 
with the fact that we are working with the trace per unit volume of an infinite, 
periodic system. 

The proof of Theorem [5] will be given in Section [S] As the proof shows, the 
theorem holds for a larger class of functions /, satisfying appropriate smoothness 
and decay assumptions. 

Our second semiclassical estimate concerns the upper off-diagonal term of 
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which we shall denote by Q!a- We shall be interested in its norm. In general, 
we define the norm of a periodic operator O by 

||0||^i = Tr [O^ (1 - /i^V^) O] . (2.14) 

In other words, ||0||^i = ||0||| + /i^||VO|||. Note that this definition is not sym- 
metric, i.e., |10||_f/i ^ ||0^|j//i in general. 
Given t as above, let us define the function 

^{p)^^go{l3{p'-fi))t{p). (2.15) 

THEOREM 3. With denoting the upper off-diagonal entry of Fa in 112. 13\} . 

we have 

\\aA - I mxM^ihV) + ip{~ihVWx))\\^, < Ch''-^^^\\i,\\H2^c) ■ (2.16) 
The proof of this theorem will be given in Section HI 



3. Preliminaries 

In this section we collect a few useful facts about the trace per unit volume. In 
particular, we recall the general form of Holder's inequality and Klein's inequality, 
both of which will be used several times in the proofs below. 

Let A and B be bounded periodic operators on either L^(R'') or L^(R''; C^), i.e., 
operators that commute with translations by a unit length in any of the d coordinate 
directions. The trace per unit volume of A is simply defined as the trace of xAXi 
where x is the characteristic functions of a unit cube, i.e., the projection onto 
functions supported in this cube. Obviously, the location of the cube is irrelevant. 
For p> 1 we also denote the p-norm of A by 

\\A\\,^{t.[A^aY^')''\ (3.1) 

Here and in the remainder of this paper, Tr denotes the trace per unit volume. We 
also use the notation ||A|joo for the standard operator norm. 
For any 1 < p < oo, the triangle inequality 

\\A + B\\^<\\A\\p + \\B\\^ (3.2) 

holds. For 1/r + 1/s = 1/p, 1 < r, s,p < oo, we have the general Holder inequality 

\\AB\\p<\\A\\r\\B\\,. (3.3) 

Moreover, if / : M x K K_|_ is a non-negative function of the form /(x, y) = 
'^igi{x)hi{y), then Klein's inequality 

Tr/(A,B) -.^Y.^^ 9^{A)h,{B) > (3.4) 

i 

holds for self-adjoint A and B. 

Inequalities p.2|) - p.4p are well-known in the case of standard traces, see [22] 
and |23]. They extend to the periodic case via the Floquet decomposition [19j 
Sect. XIII. 16]. Specifically, 
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with operating on L'^{C), and = denotes unitary equivalence. The trace per unit 
volume equals 

TrA=/ T,^. A^J^ (3.6) 

The inequalities p.2p - p.4p then easily follow from the standard ones using that 

(AB)^ ^ (3.7) 

and that g{A)^ = g{A^) for appropriate functions g and self-adjoint A. 
One also checks that 

\TtA\<\\A\W. (3.8) 
By induction, if follows from p.3p that 



\TtAiA2---A„\ < PiA2---A„||i< (3.9) 

i=l 

where 1 < Pi < oo and X^ j^r^ = ^■ 

Note that while the local trace norms introduced here share many properties with 
the usual Schatten norms, they are not monotone decreasing in p. For instance, 
if A^ is a rank one operator (whose norm is not independent of ^) then the norm 
\\A\\p is actually increasing in p. In particular, the finiteness of the p-norm does 
not, in general, imply finiteness of the g-norm for q > p. 

In the proofs below we need one more inequality that generalizes an inequality by 
Lieb and Thirring to the case of the trace per unit volume (see, e.g., [18l Sect. 4.5]. 
If A and B are periodic operators and p > 2, then 

\\AB\\p^Tt\AB\p^Tt (iB^IIAns^l)^/^ < TrlB^lf/^IAIPlB^I?'/^ (3.10) 

Again, the proof follows from the standard one using the Floquet decomposition. 
In the special case that A is a multiplication by a periodic function a{x) and B is 
multiplication in Fourier space, i.e., B = &(— iV), (I3.10p reads 

||a(x)6(-»V)||p < (2^)"''/^'||a|Up(c)||&||LP(K.) . (3.11) 

4. Properties of ao 

Recall the definition of Kt in (|1.6|) . and recall that ao denotes the eigenfunction 
of Kt^ + V corresponding to the eigenvalue zero, which we assume to be simple. 
We assume that > 0, and that V € LP{R'^), with p = 1 for d = 1, p > 1 for 
d ~ 2 and p 3/2 for d — Z. There is no parameter h in the definition of Kt] in 
other words, — 1 in this section. 

Proposition 1. 



for any k < Kf. '■= Im ^/JT+ItFT^. 

In three dimensions, this proposition also holds for n = Kc, in fact. 

The key to prove Proposition[l]is to estimate the integral kernel of K^^. It turns 
out to have the same behavior as the one for (— + k^)~^, both at the origin and 
at infinity. Hence we find it convenient to state our bound in terms of G(x — y; A), 
the Green's function of -V^ - A for A G C \ M+. Note that G{x, A) > for A < 0. 
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Lemma 1. For some constant C^^jt depending on n/T , the integral kernel of K^^ 
satisfies 

\K^\x -y)\< C^itG{x - y; -(Im + ii^Tf) . (4.2) 
Proof. We start with the series representation 

tanh^: _ 2 



which imphes that 



K-^ ^ V — (A A) 

^ ;il(n-i)2(27rr)2 + (-v2-M)2- 



We can rewrite this as 



Hence 



9 °° 1 1 

^T' = -y^Tlm n . (4.5) 



KZ^ix-y) = - ^ —^lmG{x-y;iJi + i(n-^)2TTT). (4.6) 
TT n - 5 

The result now follows easily using the explicit behavior of the Green's function G 
at and at infinity. □ 

Lemma 2. Let L he an integral operator with integral kernel bounded as 

\L{x,y)\<Ae--\--y\G{x-y;-e) (4.7) 

for K > and e>0 (or e > in d = 3). Then there is a finite constant Cp^d such 
that for any Ui,U2 e L^p{R'^) 

e'^l^lf/iXC/ae-'^l^l < ACp^de'^^^^P^-^WUilUpWahp (4.8) 

where p> 1 in d= 1, p > 1 in d=2 and p > 3/2 in d = 3. 

Proof For ip e L'^{R'^), 

< / my)\e-''\^^my)\\L{x,y)\\U^{x)\'e'^\^\ 

X \L{x,y')\\U2{y')\e-''^^'^\^{y')\dydy'dx 
<A^ [ \i,{y)\\U2{y)\G{x - y; -e)\Ur{x)\^G{x - y'; -e)\U2{y')my')\ dydy' dx 

= A^\\\U,\{-V^+e)-'\U2m\\l , (4.9) 

where we used the triangle inequality 2|a;| — lyj — |y'| < |a; — + — The result 
now follows from the fact that the operator norm of |C/i |(— V^ + e)~''^|C/2| is bounded 
by a constant times e''/^^P)^^||[/i||2p||?72||2p, as can be seen by an application of the 
Hausdorff- Young inequality (for e > 0) or the Hardy-Littlewood-Sobolev inequality 
(for e = and d = 3) and Holder's inequality. □ 
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Proof of Proposition^^ The function (f> ~ \/W\'^o is in L'^iW^) (because of the 
relative form-boundedness of V) and satisfies 

^^-VW\^^^. (4.10) 

where ^/V ^/ V I ^ I ■ For given i? > 0, we decompose (j) — 0i + 02, where 
'/'2 = <PX{\x\>R}- Introducing Ui = ~X{\x\>r}^/W\ and U2 = X{\x\>r}VV, we 
have 



i2 = f/l^f/202+/, 



where 



(4.11) 



(4.12) 



We shall now use that the Green's function of the Laplacian has the property 

G{x - y- -k") < Q/«e'(^-^)l^-^lG(x - y; -e') (4.13) 

for < e < K. In combination with Lemma [T] and Lemma [21 this implies that for 
K < Kc, the operator 

(4.14) 



is bounded. Since e^l'^'l^il < e^^l^i], we conclude that e'*'^'/ S L^(R''). Similarly, 
we observe that 



< c\\yx{\x\>R}\\t 



(4.15) 



with p as in Assumption[TJ The latter expression is less than one for R large enough. 
Hence 



e ' '02 — 



1 



= ( l-e''l^lj7i — C/2e-'^l^l 



/ 



(4.16) 



is an element of L'^i'R''^). Since obviously also e'^'^'^i G L^(IR''), this concludes the 
proof. □ 

We shall use the result of Proposition [T] in the following way. Recall that t was 
defined in (jl.l7l) as twice the Fourier transform of Kx^ctQ = —VaQ. The following 
proposition collects all the regularity properties of t and ag that we shall need 
below. 

Proposition 2. The function t in j7| ), together with all its derivatives, is a 
function in L'^{M.'^), with q — 2p/{p — 1), i.e., q — 00 for d — l,2<q<oo for 
d — 2 and q — 6 for d ~ 3. We also have that 

\dn{p)\^ 



d 1 +p2 



dp = A{x''^/\V\aQ 



1 



(4.17) 



and that J^^ (|a;'>'Vao(x)p + \x'^ao{x)\^) dx < 00 for all 7 G Nq. 



Proof. The L'^ property follows easily from the Hausdorff- Young [17, Thm. 5.7] and 
Holder inequality, 

l|a^i||2p/(p-i) < 2\\x^Vaohp/ip+i) < 2\\V\\l/^x''^\ao\\2 , (4.18) 
using Proposition [1] and our assumption that V ^ L^. 
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The boundedness of (|4.17p follows again from Proposition [T] and the fact that 
the operator on the right side is bounded by our assumption on V. Finally, the 
last statement follows easily from (j4.17p writing the integral in Fourier space and 
using that ao{p) = t{p)/{2KT^{p)), where Kx^p) denotes the function obtained by 
replacing — by in (|1.6p . □ 

For d = 3, one could even drop the 1 in 1 — and 1 in (|4.17p : boundedness 
of the operator follows from the Hardy-Littlewood-Sobolev inequality. 



5. Proof of Theorem [T] Upper Bound 

In this section, we shall prove the upper bound in Theorem [1] i.e., we shall show 
that e < Ch with e defined in ([1121) . We shah denote /3 = l/T and = l/Tc. 
RecaU that P = I - h"^ D with D > 0. 

By definition, F^^^{T,n) < T^'^^(T)~T^'^^{Ta) for any admissible state P. As 
a trial state, we use 



Ta^I 2^ 1= ^ (5.1) 

where Ha is given in (|2.5[) with A as in (j2.1l) . Note that Pg, defined in (|1.9|) . indeed 
corresponds to setting A = in ()5.1|) . For t, we choose (|1.17p . which is refiection 
symmetric and can be taken to be real. We normalize it such that (|1.18|) holds. 
The integral kernel of A is then given by 

M^, y) = (IJ^ W^) + Hy)) V{h~\x - y))ao{h-\x - y)) . (5.2) 

Note that Ha is unitarily equivalent to —Ha, 

UHaW^ = -Ha with U= J ^ . (5.3) 

Hence also UTaU^ = 1 — Pa and, in particular, 

SiTA) - -iTr [Pa IuPa + (1 - Pa) ln(l - Pa)] ■ (5.4) 

A simple calculation shows that 

Pa In Pa + (1 - Pa) ln(l - Pa) - TolnPo - (1 - Po) ln(l - Po) 
= -PHaTa + PHoTo - In (1 + e'^^^) + In (l + e-^"°) . (5.5) 



1 



Moreover, 



HATA-Horo=[ fc(7A-7o) + AdA fcaA + A(l - ^a) \ (^ g^ 
' A7A + kaA k{lA ~ 7oj + A^a / 



MICROSCOPIC DERIVATION OF GINZBURG-LANDAU THEORY 



15 



where k denotes the left upper entry of Ha (and Hq). From ()1.4[) and (|5.4p - ()5.6p 
we conclude that 



_ 1 



Tr [ln(l + e-^^^) - ln(l + e-^"°)] 



Cxi 



2 dxdy 
(2^ 



+ 



CxI 



dxdy (5.7) 



where the Tr in the first term on the right side has to be understood as the sum 
of the traces per unit volume of the diagonal entries of the 2x2 matrix- valued 
operator. In general, the operator ln(l + e^'^^^) — ln(l + e^'^^°) is not trace class if 
A is not, as can be seen from (p3]) - ((?^ . In the evaluation of J"^^^(rA)-^^^^(ro) 
only the diagonal terms of (|5.6p enter, however. 

The first term on the right side of (|5.7I) was calculated in Theorem [2] above. 
Since /3 = /3c + 0{h?), we can replace /3 by /3c in all the terms of order /i^, yielding 
an error of order . For the term of order h? , we obtain 

h^P I t{qfgo{P{q'~y))dq 



JR'' 

+ h^PcD [ t{qf (5o(/3c(9' - m)) - Pc{q^ - M).9i(/5c(g' - m))) dq + 0{h') 



where D = {Tc — T)/{h'^Tc). The term in the last line of order equals 



D 



tiqycosh-'i^l3,{q'-fi))dq = 



3i(/3c(g'-/i)) 



9^ - M 



(5.8) 



dq (5.9) 



according to our normalization (|1.18l) . 

The second term on the right side of (|5.7p can be rewritten as 



Rt^xC 



Vih-\x - y)) \'^{^{x) + i,{y))a^{h-\x - y))\ dxdy 



pe(27r 



16 



x)\ao{x)\ cos (|p • x) dx 



t{q) go(/3c(g' - m)) (2i(g) + t{q - hp) + t{q + hp)) dq . 



By writing 

2t{q) + t{q - hp) + t{q + hp) 

= At{q) + h" [{p ■ ^)h] (q) + 



(5.10) 



[{p ■ V)H] {q + shp){l ~ |s|)3 ds (5.11) 
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we observe that ()5.10|) equals 



tiqy goWciq- fJ-)) dq 



E i£>^^\9,i^) I til) [^^^,t]{q) goWdq^ - m)) dq + Oih^'") , (5.12) 
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where the error terra is bounded by 



Ch' 



6-0 



if2 



\V{x)\\ao{x)\^\x\Ux. 



(5.13) 



The latter integral was shown to be finite in Proposition [TJ 

If V < 0, the last term in (15. 7p can be dropped for an upper bound, but we 
do not need to make this assumption. Since V is relatively bounded with respect 
to the Laplacian, we can bound the term by an appropriate norm. Recall the 
definition of the norm of a periodic operator in (|2.14p . For general periodic 
operators O, we have the bound 



V{h-'{x-y))\Oix,y)\^ dx dy 



CxI 



< 



(1 - h^\/^)-^^^V{h-\ ■ ~ y)){l - /i2v2)-i/2 ||0||2^^ 

(i-v2)-i/V(-)(i-v2)-i/i \\o\\j,^. 

The operator in question can be written as 

O = QfA - I {ip{x)So{-ihV) + So(-i/iV)'0(a;)) . 



(5.14) 
(5.15) 



Recall the definition of (p in (|2.15l) . It equals q?o for j3 ~ Pc- Hence 

O = ttA - I {ip{x)ip{-ihV) + ip{-ihV)ip{x)) 

+ J {4'ix)r){~ihV) + ri{-ihV)ij{x)) , (5.16) 

with 

viq) = {Ml3{q^ - m)) - PcgoiMq^ - m))) t{q) . (5.i7) 

For the term in the first line of (|5.16p . we can apply Theorem [3] to bound its 
norm by C/i^~^/^||'0||_ff2 ■ The norm of the term in the second line of (15. 16^ can 
be bounded by 

1/2 



(5.18) 



It is easy to see that \ri{q)\ < C(/3 - l3c)\t{q)\/{l + q^), hence ((^?T51) is bounded by 

C/l3-'i/2||y,||^i. 

For ■(/;, we shall take a minimizer of the GL functional (|1.7I) . Under our Assump- 
tion [5] on W and A, it is easily seen to be in . For this choice of V', we thus 
have 



F 



BCS/ 



(r,(i - h^D), ^) < -F^^S(rA) - i^^^^(ro) < h^-"^ {e^^ -B3 + ch) (5.19) 



for small h. This completes the proof of the upper bound. 
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6. Proof of Theorem [T] Lower Bound, Part A 

Our proof of the lower bound on ^-j Theorem [1] consists of two main 

parts. The goal of this first part is to show the following. Let again Pq denote the 
normal state defined in (|1.9p , which is the minimizer of J^^^^ in the non- interacting 
case V = 0. We claim that for any state P satisfying J'^^^(P) < J'^^^(Po), we can 
decompose its off-diagonal part a as 

a= |(V'(x)So(-iW) + ao(-iW)V'(.T)) +^ (6.1) 

for some periodic function ip with H^{C) norm bounded independent of and with 
UWh^ < 0(p-'^''^), where we use again the definition (I2.14[) for the norm of a 
periodic operator. This latter bound has to be compared with the norm of the 
first part of which is 0{h^-'^/^) (for fixed V ^ 0.) 

The remainder of this section contains the proof of (|6.ip . It is divided into three 
steps. 

6.1. Step 1. We claim that for any state P of the form (|1.2p satisfying T^'^^{T) < 
jrBCS(p^)^ we have that 

^ Tr (aa)2 + ^(«( • , y)\K^-'^ + V(h-\- - y))\a{ ■,y))dy<0. (6.2) 

Here, K:^'^ denotes the operator 

A.w i-ihV + hA{x)f -^i + h^W{x) 

= /a /. ~ .W . (D-^J 



tanh 



f ((-iW + hA{x)f - M + h^Wix) 



with /3 — 1/T. In (|6.2p . it acts on the x variable of a{x, y), and ( • | • ) denotes the 
standard inner product on L^(E'^). Note that k!^'^ differs from the operator Kt 
defined in (|1.6p by a scaling by h. For T — Tc, the ground state of k!^'^ + V{h^^{ ■ — 
y) equals h^'^^^ao{h^^{x — y)) (up to an /i-independent normalization). 

Using that {—ihV + hA{x)) — ^ + h'^W{x) = Tln((l — 7o)/7o), we may write, 
for any state P, 

_^BCS(r)„_^BCS(r^)^ iy^(r,Po)+ / V{h-\x-y))\a{x,y)\' dxdy , (6.4) 

where % denotes the relative entropy 

H(P,Po) = Tr [P(lnP-lnPo) + (l-P)(ln(l-P)-ln(l-Po))] . (6.5) 
We have the following lower bound. 

Lemma 3. For any < P < 1 and any Po of the form Po = (1 + e")-'^, 



H(P,Po)>Tr 



(r - Po)^' 



^Tr [P(l-P)-Po(l-Po)]' . (6.6) 



tanh(i?/2) 

A similar bound as (|6.6p . without the last positive term, was used in [10] . 
Proof. It is elementary (but tedious) to show that for real numbers < < 1, 

In ^ — ^ 

x\n- + {l-x)\n\^ > ^^(x-j/)2 + l(x(l-a;)-y(l-y))' . (6.7) 
y 1-2/1- 2y 3 

The result then follows from Klein's inequality (|3.4p . □ 
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Note that for our Fq, H equals /3i?o, which is diagonal as an operator- valued 
2x2 matrix. Hence also Hq/ tanh(/3iIo/2) is diagonal. Its diagonal entries are 



j3 K^'^ and (3 K^'" , respectively, where K^'^ is given in (16.31) above. Hence 



^-rTA,W 



Tr 



tanh i^Ho 



(r-To) 



2Tr 



(7 - 70)^ 



2Tr 



T^A.W - 



Since x/ tanh(a::/2) > 2, we can replace K:^'^ by 2T for a lower bound. We shall 
use this in the first term on the right side of (|6.8p . 
For the last term in (16. 6|) . we use 

(6.9) 



Tr [F(l-F)-Fo(l-Fo)]' >2Tr - j) - ^^{1 - j^) - aaf 



We claim that 

2 Tr (7 - 7o)^ ['^'^-^ ~ 7) - 7o(l - 7o) - aaf > ^ Tr {aaf . 

This follows easily from the triangle inequality 

||aa||2 < 117(1 - 7) - 7o(l - 7o) ^ ""lb + ||7(1 - 7) - 7o(l - 7o)ll 
together with the fact that 

Il7(l - 7) - 7o(l - 7o)ll2 < II7 - 70II2 , 

which can be seen using Klein's inequality p.4p . for instance. 
This completes the proof of 



(6.10) 



3.11) 



3.12) 



6.2. Step 2. Recall that K!^''^+V{h ^( • )) is non-negative and has a non-degenerate 



isolated eigenvalue zero. Hence it will be convenient to replace K:^'^ by K^'" in 
(|6.2p . The following lemma quantifies the effect of such a replacement. 

Lemma 4. For T = Tc — 0{h^) and h small enough, 



K, 



A,W 



Vih-\- -y))>^(^K'^^ + V{h-\- -y))) 



Ch^ 



(6.13) 



for a constant C > depending only on ||VF||oo, ll^llci '^'^'^ h ^(T— Tc). 



The proof shows that the prefactor 1/8 can be replaced by any number less than 
one, at the expense of an increase in the constant C. 



Proof. We start by noting that [T, (4.3.91)] 



tanh(a;/2) 



^ a:2/4 + 



fc27r2 



E 

k=l 



2fc^7r 



2.^2 



(6.14) 



In particular, this function is operator monotone as a function of x^. Let p = 
—ih\7 + hA{y) and A{x) — A{x) — A{y). Since the ground state of the operator 
Rrp'^ -j- V{h~^{ ■ — y)) is localized within a distance 0{h) of ?/, A{x) is, as far as the 
ground state is concerned, effectively a perturbation of order h since A is assumed 
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to be Lipschitz continuous. Using Schwarz's inequality, 



(p + hA{x)y - fi + h^Wix) 



>(l-e) + hp- A{x) + hA{x) ■ p- fi 



2 



\W + A 



2||2 

oo 



>(l^er[p'-^,]'--\\W + Am--\p-Aix)+A{x)-py 
We write the right side as (3^^{R~ Q), where 



Q 



p ■ A{x) + A{x) ■ p 



and 



,2/1 ,^2 M ^''/^^llw 1 72||2 



R^fi'il-erip'-f,]' 



\W + A^ 



3.15) 



3.16) 



(6.17) 



As long as i? — Q + Att'^ > 0, we can use the operator monotonicity of (|6.14p to 
obtain a lower bound. This condition is certainly satisfied for small enough values 
of e, h'^/e and /i^/e. Hence we can use the resolvent identity to conclude that 



oo 



R 



2 ^ /c27r2 + i?/4 

fc2^2 ^ 1 



2 k^TT^ + R/4^k^n^ + R/A 

S 1.2^2 I R/zl'^i 



8 A:27r2 + i?/4^ fc27r2 + {R-Q)/A^ k^ir^ + R/A 

Since also R — Q > Q — Att"^ for e, /i^/e and /i*/e small enough, we can bound 

1 



3.18) 



g- 



P7r2 + (i?- 0)/4 



Q<AQ 



for all fc > 1. We thus obtain the lower bound 

R 



1 / 1 ^ 

^ V ^ 2 ti '^^'^ 



^+R/Al 13 fc27r2 + i?/4 ^ fc27r2 j^jij^ 



(6.19) 



(6.20) 



We start with deriving a lower bound on the first term on the right side of (|6.20p , 
and defer the discussion of the second term to (|6.28p et seq. For h'^/e small enough, 
the first term on the right side of (|6.20p is bounded from below by 



(1 - efl3K, 



e 12 7r2-/i4e-i/32||^ + ^2|j2^ ' 



-(3^\\W + A'\ 



3.21) 



where we denote by Ay the constant vector potential A{y). Moreover, it is elemen- 
tary to show that 



> K^"' - 2{Tc - T) 



3.22) 
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for T <Tc. We further have 

(1 - 26)<;'° + V{h-\- - y)) = (1 ~ 36) + V{h-\- - y))) 

+ e(<-° + 3n/.-i(. ~y))) 

>(l-3e) {Kp/ + V{h-\- -y)))-Ce, 

(6.23) 

A 

where — C is the ground state energy of Kj,^' + iV, which is bounded by our 
assumptions on V . It is also independent of Ay since Ay may be replaced by zero 
by a unitary (gauge) transformation. 

We now want to get rid of the constant vector potential A{y) in Kj,^''^ on the 
right side of (|6.23p . We claim that 

+ V{h'\ ■-y))>l (<f + Vih-\- - y))) - Ch'Aiyf . (6.24) 

This follows from the fact that 

> ^K^T° + IkI^''° - Ch^A{yf (6.25) 

(since the function (j? — tanh[|-(p^ — /i)] has a bounded Hessian) and 

K%^J'^ ^V{h-^{- -y))>^ (6.26) 

by the definition of T^. 

For 3e < i and Tc - T = 0(/i^), we conclude that 

K^^^ + V(h-\- ^V))>\ {K'^f + V{h-\- - y))) 

-C{h^ + e + h^e-^)-E (6.27) 

where E denotes the last term in (I6.20p . Let P denote the projection onto the 
ground state of X^'^ + V{h~^{ ■ — y)), given by the function ao{h~^{ ■ — y)). Let 
also P'^ = 1 — P. Since E is positive, it follows from Schwarz's inequality that 

E < 2PEP + 2P'EP'' . (6.28) 

From the assumption on A G it follows easily that < Q < Ch^e^^{l +p^). 
This immediately implies that E < Ch'^e-'^{1 + jP) < Ch'^e-^{1 - h^V^) (since A 
is bounded). We shall choose e = O(h^), with h'^/e small enough to ensure that 

i (^K^^ + V{h-\ ■ - y))) - 2P''EP'' > 0. (6.29) 

That this can be done follows from the above bound on E and the fact that k!^'^ + 
V{h-^{- -y)> vP'^il - h'^V^)P'' for some J/ > 0. 

It remains to show that PEP < 0{h^). Since P has rank one, we can write 

PEP = Ch-'' {aoih-\ - - y) \E\ aoih-\ - - y)) P (6.30) 



MICROSCOPIC DERIVATION OF GINZBURG-LANDAU THEORY 



21 



for some /i-independent constant C determined by the normalization of ao- We 
have 

Q = ^^-^ {2p-A{x)+ih<:\hjA{x)^ {2A{x) ■ p - ih <1\\ A{x)^ 

< ^ (8p.I(a;)I(a;) •p + 2/i2(div^(a;))') . (6.31) 

Note that the last term is bounded by our assumption A G C-^. From the Lipschitz 
continuity of A it follows that |A(a;)| < C\x — y\, and hence 

p- A{x)A{x) -p^C^p-lx-yl^p. (6.32) 

We thus have PEP < Ch'^Pj^k ^^^fc '^here 



i—1 



p, + hAi{y) ^ 

' V \ 9 "0 \P) 

Pp7r2 + 1/32(1 _ e)2 [(p + /iAy)2 - _ ^ 



dp 



with 8 = 0{h^) denoting 6 ~ |;/i^/3e^^||M^ + For small enough h, we can thus 
bound 

Xk < C— / (1 + \xf) (|Vao(a;)|' + \ao{x)\^) dx . (6.34) 
eft jRd V / 

The latter integral is finite, as proved in Proposition[2]in Section^ Hence k^Xk 
is bounded, uniformly in h for small h. This completes the proof. □ 

6.3. Step 3. In combination with Lemma |4l we conclude from (|6.2p that for T = 

4T 1 r 

— Tr[aap + - {a{- ,y)\K'^^ + V{h'\- - y))\a{- ,y)) dy < Ch^Ml (6.35) 



for any state F with J^^^^(F) < J"^^^(Fo). We shall now show that this inequality 
implies (|6.ip . 

Recall that the operator if^'^ + V[h~^{ ■ — y)) on L-^iW^) has a unique ground 
state, proportional to aQ{h^^{x — y)), with ground state energy zero, and a gap 
above. Normalize ao as in (|1.18p . and let 

i^{y) = {2tiY'^ (h ( \ao{x)\'' dx\ ( ao{h-^[x-y))a{x,y)dx. (6.36) 
Note that i/i is a periodic function. If we write 

oi{x,y) =il){y)-^—yj:^ao{h ^[x ~ y)) + ^o{x,y) (6.37) 

the gap in the spectrum of + V{h~^{ ■ — y)) above zero, together with (|6.35l) . 
yields ||^o||2 < 0(/i)||Q;|j2. We can also symmetrize and write 

a{x,y) = l{-^{x)+^{y))^-^ao{h-\x~y))+^{x,y), (6.38) 

again with ||^||2 < 0(ft.)||a||2. In order to complete the proof of (|6.ip . we need to 
show that llV'lli/i is bounded independently of ft., and that the norm of ^ is 
bounded by C'(/i2-'^/2). 
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An application of Schwarz's inequality yields 



/fid \ao{x)\'^dx 



Moreover, 



lall^ < 



'2 ^ {2^)d 
Since |l^o||2 < 0(/i)||q:||2, we see that 



\4,{x)\''dx / \ao{x)\''dx+\M 



•2-d 



\il'{x)\'^dx I \aQ{x)\^dx. 



12 ^ V- ' -V- ,,^^^Y 

Again by using Schwarz's inequality, 

\T7,( M2^ ^10 ^rf.rf-2 /Ilt^xcK^^+^^/)"(^'y)l' ^^^^^ 

|V^(.)| dx < (2.) h J\a,ixWdx • 

In order to bound the latter expression, we use the following lemma. 
Lemma 5. For some constant C > 0, 



(6.39) 
(6.40) 
(6.41) 

(6.42) 



'xC 



\{Vx + Vy)a{x,y)\ dxdy 

<C j^{a{ - ,y)\K°^^ + V{h-\ - - y))\a{ - ,y)} dy (6.43) 

for all periodic and symmetric a (i.e., a{x,y) — a{y,x)). 
Proof. By expanding a{x,y) in a Fourier series 

a{x,y)= e^P-(^+?^)/25p(a; - y) (6.44) 

and using that 5p(x) — ap{—x) for all p E (27rZ)'^ we see that (|6.43p is equivalent 
to 

K}^'° + K^}"-" + 2 Vix/h) > ^h^p^ (6.45) 

for all p e (27rZ)''. This inequality holds for all p G M'', in fact, for an appropriate 
choice of C > 0, as we shall now show. 

Since K^^'^ > C(l + h'^{—iV +p/2)^), it suffices to consider the case of hp 
small. Let k denote the gap in the spectrum of k!^° + V{h^^ ■ ) above zero, and 
(Pq{x) := h~'^^'^(po{x/h) its normalized ground state. (Note that ipo is normalized 
to one, and hence equals a constant times ag.) Then 



K, 



ip,0 



K, 



-ip,0 



> K 

> K 



2V{x/h) 



^•P/2(l-|^5)(^S|)e 



— ix-pl2 I —ix-pjl 



1 - 



\^Q{x)\^e-'^''-^dx 



(6.46) 



where the last expression is simply the lowest eigenvalue of the operator on the 
previous line. Since (^o is reflection symmetric. 



/ 



|(po(a;)Pe-*''^-Pda; 



|i^o(a^)| <^os (/ix • p) da; = 1 — 0{h p ) 



3.47) 
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for small hp, since J^^ \x\'^\(po{x)\'^dx is finite by Proposition [2] 



□ 



By combining (|05)) with (jOTI) . ((O^ and we see that || VV'lb is bounded 

by a constant times HV'lb- To conclude the uniform upper bound on the norm 
of ip, it thus suffices to give a bound on the norm. To do this, we have to utilize 
the first term on the left side of Eq. (j6.35|) . 

Recall that a can be decomposed as a = hao^) + £.0, as in ()6.37|) . where ao is 
short for the operator ao{~ihV). The following lemma gives a lower bound on 
(Tr (aa)^)^/'*, the 4-norm of a. 



Lemma 6. For some < C < 00 we have 



|a||4> 



\tl;{x)\'^dx / ao{q) 



dq 
(2^ 



ch'-'M% 



-,1/4 



1(C) 



where 



^ Ch^-''/^\\i;\\l^^ (l + Ch^-^/^yj\\i 
max{0, • } denotes the positive part. 



1/2 



(6.48) 



Proof. By the triangle inequality, 

II^IU > h\\aa^p\\4 - ll^olU . 



(6.49) 



We can bound the last term as W^oWl < ||?o||oo||^o||2- Recall that we have already 
shown that H^olb < 0{h)\\a\\2 < 0{h'^-'^^^)\\i}\\2. Note that we cannot bound the 
00-norm simply by the 2-norm, since the norms || • ||p defined via the trace per unit 
volume in (13.11) are, in general, not monotone decreasing in p. We can, however, 
use that ||Co||oo < \\cf\\oo + /iUccoV'lloo and that ||a||oo < 1 for any admissible state. 
Moreover, we claim that 

||«o^lU < C,U\Uh-^/^ (^J^^ \aoix){l + |a;|)n'/' dx^ ' (6.50) 

for V > d\ that the latter integral is finite follows from the last statement in Propo- 
sition[21 which says that ao( • )(1 + | • is in L^(E'*) for any v > ^. 

Eq. (|6.50p can be obtained with the aid of Young's inequality [Tjl Thm. 4.2], as 
we now explain. With x denoting the characteristic function of the unit cube, and 
Xj(a;) = x(a; — 3) for j € , we have, for any /, g € L^(M'^), 



l(/l«o^l5)l 

< E (1 



[b-fc|-v/rf]H 

h 



X 



|x,(x)/(x)ao(^)(l 



h 



^{y)xk{y)9{y) 



dx dy 



<ch-'/^ J2 (1 



■\]-k\-Vd]+ 

h 



IXj/ll2||Xfc5ll2||V'llL'i(c)l|ao(l 



irii4/3, 

(6.51) 



where we used Young's inequality in the second step. We can further bound 
IIXi/l|2||Xfeff||2 < A||xi/||2 + A^lxfeslli for A > 0. After doing the sum over j 
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and k and optimizing over A, we thus have 

\{fM\9)\ 



<Ch-'^/nfh\\g\\2j2[^ 

jei,'' 

which yields (|6.50p . 

We have thus shown that 



h 



ao(l 



l4/3- 



3.52) 



ll^oll^ < Ch^-'/^^A\2 (l + Ch'-^^^Mi) . (6.53) 
It remains to investigate the first term in (I6.49p . A short calculation shows that 
\\ao^A\t = h-'' i^iPiWiP2MP3)r{-pi-P2-P3)Fihpi,hp2,hp3), (6.54) 

Pl,P2,P3 

where 

F{pi,p2,P3) ^ / ao{q)aoiq+Pi)aoiq+Pi+P2)ao{q+pi+P2+P3) -77m ■ (6-55) 
Note that F is bounded. In fact, 

\Fipi,p2,P3)\ < (2^)-''||ao|it (6.56) 

To see that this is finite, one can use that |So(p)| < C\t{p)\/{1 + p"^) (from the 
definition (jl.l7p ) and that t ^ L'^ with q as in Proposition [2] In the same way, one 
can show that F has a bounded derivative, and hence 



\F{pi,P2,P3) - i^(0, 0,0)1 < C(|pi| + \P2\ + \P3\ 

We are left with giving a bound on 

^ H{piW {p2)'<P{P3)i'* {-Pl ~ P2 - P3) \Pl\ 



(6.57) 



Pl,P2,P3 



< 



Y.\^{p)?\p?] E 



\Vil;{x)\'^ dx 



1/2 



E yi'* {P2)'4'{P3W {-P - P2 - P3) 

P2,P3 

1/2 

\'ipix)fdx' 



1/2 



(6.58) 



where ip is the function whose Fourier transform equals {ipl- For d < 3, we can use 
the Sobolev inequality 

\ 1/3 

mx)\'dx) <c\mi^c) = c\mmic)- 



(6.59) 



Hence 



dq 



\\ao^\\t>h-'' mx^dx / aoiqr-^~Ch'-'Urm^c)- 
Jc JR<i i^^r 

and this completes the proof. 



(6.60) 

□ 



We have already shown that ||V?/'||2 < C||?/'||2, which also implies that HV'IU < 
CllV'lh via Sobolev's inequality for functions on the torus. If we use also that 
IIV'IU ^ llV'lh (since the norms of ip are defined via integration over the unit cube 
C), we conclude from (jOSll that ||a||4 > Ch^-'^/^iWtPh - CUWl^^) for h small 
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enough. In combination with (|6.35p and (|6.4ip this implies that \\ip\\2 < C This 
shows that the norm of ijj is indeed uniformly bounded. 

It follows that 11^112 < 0(/i2~'^/2). To conclude the proof of (|5T^ . we need to 
show that also ||^||hi < 0{h^-'^^^), i.e., ||VC||2 < 0{h^-'^/^). We have 

^{x,y) ^ Ux,v) + \ {^{x) - i;{y)) j^-^ao{h-\x - y)) . (6.61) 

From the definition (|6.37p it follows easily that ||Co||_ffi < 0{h'^^'^/^), since 

Uofm <C j^{U-.v) |^T° + V{h-\ - -y)) + l\U-. v)) dy 

<C{h^a\\l + \U\D<0{h'-'') (6.62) 

using (j6.35p . Recall that the definition of the norm in (|2.14|) is not symmetric, 
hence this does not immediately imply a bound on the norm of ^. However, we 
can estimate 

/" |VV'(x)|2 / \ao{h-\x-y))\^dxdy<0{h'-'^). (6.63) 
Jc Jm.'i 

Finally 

h-^" [ \i;{x)-i,{y)\'\\/ao{h~\x-y))\^dxdy 

= 4/1-'* J2 IV'(P)P / |Vao(a;)p sin^ (i/ip ■ x) dx < 0{h'^-'^) , (6.64) 
pe(27i-z)<i 

since the norm of ^ is bounded and J jVaoPl^^piia; is finite, as shown in Propo- 
sition [2] in Section |4l This proves the claim. 

7. Proof of Theorem [T] Lower Bound, Part B 

We now conclude the proof of the lower bound of Theorem [TJ For convenience, 
we shall divide the proof into 2 steps. 

7.1. Step 1. Let F be a state with F^^^{T) < J^^'^^{To). In the previous section, 
we have shown that the off-diagonal part a of F can be decomposed in the form 
(|6.ip . Given ijj defined in (j6.ip and some e > 0, define ■(/;< via its Fourier transform 

$<ip) ^ me{eh-' - \p\) , (7.1) 

where 9{t) = 1 for t > 1, and otherwise. The function is thus smooth, and 
||'0<llff2 ^ C(l + £h^^) since V' is bounded in H^. We shall choose h < e < 1. 

Let also 4'> = ^ ^ Since ■0 is bounded in H^, the Lp'iC) norm of V'> is 
bounded by 0{he^^). We absorb the part i(V'>(a;) + 'il)^{'y))aQ{h~^ {x — y)) into ^, 
and write 

a{x,y)^ ^{'ilj<{x)+i;<{y))j^--^ao{h ^{x - y)) + a{x,y) (7.2) 



where 



a{x,y)=^{x,y) + ^{'ip>{x)+il;>{y))j^-:^aoih '^{x-y)). (7.3) 
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In the previous section, we have shown that HCIIffi < 0{h^ '^^^). From the bound 
||'0>||2 < 0{he-^) it thus follows that \\a\\2 < Oik'^^'^^^e^^). Note that we cannot 
conclude the same bound for the norm of a, however. 

As in (|2?T|) . let A = -i(V'<(a;)^(-iW) + t{-ihV)ip^{x)). Its integral kernel is 
given in (15.21) . with -0 replaced by Let H/^ be the corresponding Hamiltonian 
defined in (12.51) . We can write 



^ -Ltt [ln(l + e-^"^) - ln(l + e-l^"")] 

- h'-'" I V{h~\x - y))\ \^^^{x) + V<(2/)|' ~ y))\^ ^"^^^ 



+ \T'H{T,T^)+ f V{h-\x~y))\a{x,y)\^dxdy, (7.4) 

JcxR'' 

where H denotes again the relative entropy defined in (|6.5p . 

The terms in the first two lines on the right side of (|7.4p have already been 
calculated. The first term is estimated in Theorem [21 and a bound on the second 
term was derived in Section [5] on the upper bound. The error in replacing (3 with 
/3c is as for the upper bound. Using the uniform upper bound on the norm of 
'0<, as well as ||V'<||if2 < Ce/h, we obtain the lower bound 



+ lTn{T,rA)+ / Vih-\x-y))\aix,y)\'dxdy. 

(7.5) 



It remains to show that the terms in the last line of (|7.5p are negligible, i.e., of higher 
order than h'^~'^, for an appropriate choice of e ^ 1. This will be accomplished in 
the next step. 



7.2. Step 2. We again employ Lemma |3] to get a lower bound on the relative 
entropy 'H(r,rA). It implies that 



TH{T,Ta) > Tr 



Ha 



tanh hfiHA 



(L-Fa)^ 



(7.6) 



As in the proof of Lemma IH we shall use the fact that x i— > y/x/ tanh-^x is an 
operator monotone function. 

An application of Schwarz's inequality yields 



Hl>{l- fi)Hl 



(7.7) 



for any 77 > 0. To bound we proceed as in the proof of Lemma SI specifically 
using ((05|) and ([S3T|) . which states that [-iV A{x) - iA{x) < C(l - V^). The 
choice e = 0{h) in (|6.15p yields the lower bound 



Hl> {I - 0{h))[-h'^V'^ - n\ 



-Ch. 



(7.8) 
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The operator monotonicity thus hnpHes that 

2 < ■ 



K. 



0,0 



(1 - 77 - Oih))-^/WHi+V-'mi. + Ch 



tanh 



1^(1 _ ^ _ o{h))-^/^yTiTV^ML + ch 



<{l-il-0{h)) 



-1/2. 



77-^11 AllL + C/.) 



tanh 



<{l~V- 0{h)r^l^ (1 + 1/3^ (r,-i|| A||L + Ch)) 
for < 77 < 1. With the choice 77 = 0(||A||oo) this gives 



tanh ^(iHA 



Ha 



> (l-0(/7 + ||A||oo))<'°®Ic^ 



tanh 1/3 Ha 



In particular, we infer from (I7.6P that 
T 



(7.9) 



(7.10) 



(7.11) 



:W(r, Ta) > (1 - 0{h + II Alloc)) Tr ^K'^^ia - aA){a ~ cia) 

where aA denotes again the upper off-diagonal entry of Fa. 
From the definition of A, we see that 

||A||co < /l||V'<||oo||t||oo. (7.12) 

Moreover, since the Fourier transform of ■(/;< is supported in the ball \p\ < e/h, 

1/2 



< 



j:\^<iP)\<\\^<\\HHc)( E 

P \\p\<eh-^ ^ J 



1 



for d = 1 



< C X <^ y^ln{e/h) for d = 2 
[ \/e/h for d = 3. 

Recall the decomposition (|7.2p of a, and define (f) by 

OA = I {ip<{x)ao{-ih\7) + ao{-ihV)^<{x)) + 4> . 

We thus have 

a — a A — <7 — (f) . 



(7.13) 



(7.14) 



(7.15) 



Since ||7/'<||ff2 < 0{e/h), Theorem [3] implies that ||0||jji < 0{eh'^-'^/'^). From the 



positivity of ° we conclude that 



Tr 



<^°(a-0)(5-0) 



> Tr K'^ji^aa - 2 Re Tr X'^^ac, 



(7.16) 



The terms quadratic in a are thus 

(1 -(5)TrX^'VCT+ / y(/i-i(a;-y))|cr(x,2/)|2dxd?7 (7.17) 

with (5 = 0{h + II A||oo). Pick some 5 >Q with (5 + 5 < 1/2, and write 



(1 - ^)if^'° + y = 5if 



T 

0,0 



1 - 25 - 25] { /4'° + ) + ( <5 + <5 



2i:ir,/i2-CU + (5 , 



+ 2^ 



(7.18) 
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where we have used that V is relatively form-bounded with respect to K^'^ to 
bound the last term, and > - 2{Tc - T) = - 2h'^DTc to bound the 
second. Using also that Kj,' > —h^V^ — fi, we thus conclude that 



(1 -(5)TrX^'' Vct+ / V{h-\x -y))\<j{x,y)\^ dxdy 

> S{Ml, - (C + M + l)lk||^) - {CS + 2DTch^) Ml . (7.19) 

Recall that \\a\\2 < 0{h^^'^/'^ / e). We shall choose (5 = if the first parenthesis on 
the right side of (|7.19p is less than (and, in particular, if it is negative), 

while 5 = 0(1) in the opposite case, i.e., when Hcrlllf^ > 2(C + /i + l)||cr|||. In the 
latter case we shall have the positive term (5||(t||^i/2 at our disposal, which will be 
used in (17^ below. 

We are left with estimating the last term in (|7.16p . It can be bounded by the 
product of the norms of a and 4>. This turns out not to be good enough, 
however. Recall from (|7.3p that cr is a sum of two terms, ^ and cr — f , where the 
latter is proportional to ijj^, and ||^|j_ffi < 0{h'^^'^^^). Moreover, as the proof of 
Theorem [3] in Section 12] shows, </> is the sum of two terms, rji and (f> — rji, with rji 
defined in (|9^ (with -0 replaced by -0<) and ||0 - ?7i||ffi < 0{h^~'^/'^). Now 

TrA'°'"(CT-^7;i =0 (7.20) 

as can be seen by writing out the trace in momentum space and using that ■;/'< and 
'!/'> have disjoint support. Hence 

Re Tr<'°a0 < C (||e||ifi + M HA\<t> - h^) 

< Oieh^-"^) + 0{h^-^/^)\\a\\m ■ (7.21) 

In the case ||cr||/fi < C||(t||2 (corresponding to (5 = above) we can further bound 
llo'llffi < 0{h?^'^/'^ / e). In the opposite case, where 5 = 0(1), we can use the 
positive term (5||cr||^i/2 from before and bound 

^-Ml, 0{h^-'''^)MH^ > -0{h^-'') , (7.22) 

which thus leads to an even better bound. 

In combination with (17.51) these bounds show that 



j^d-i (^BCS(p) _ ^BCS(p^)) > f °^(V'<) - S3 - Ce (7.23) 

where 



h h 



for d = 1 



e = /i + e2 + £+ - + — X { yjH(./h) for d = 2 (7.24) 

■\/e//i for = 3. 

The choice e = h^'^ for d = 1, e = h^/^[\n{l/h)Y/^ for d = 2 and e = h^^^ for d = 3, 
respectively, leads to 

r for d = 1 

e<Cx\ /ii/3[in(i//j)]i/6 for ^ = 2 (7.25) 
[ for d = 3. 

The completes the lower bound to the BCS energy. The statement (|1.15p about 
the minimizer follows immediately from (|7.23p and (|7.2p . 
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8. Proof of Theorem [2] 

For simplicity, we prove Theorem [2] only in the case d = 3. The cases d = 1 and 
d = 2 are very similar and are left to the reader. 

Note that the function / in (I2.6P is real analytic and has a bounded derivative. 
The second and higher derivatives decay exponentially. In particular, the following 
lemma applies. It will be used repeatedly in the proof of Theorem [2j 

For a general smooth function / : M — > K, let [oi, . . . , oat]/ denote the divided 
differences [3j , defined recursively via 

N/ = /(a) : [ai,a2j/ = (8.1) 

ai — 02 

and 

r „ 1 [«!,..., flAT-l]/ ~ [a2,...,aN]f f„ „^ 

is case ai ^ a^. The extension to coinciding arguments is simply by continuity. In 
case of distinct arguments, 

N 

[ai,a2, . . .,aN]f = ^ /(a^) Y[ i^j - Oi)"^ ■ (8-3) 

In case / is analytic in a neighborhood of the real axis, we have 

If ^1 
[a„...,a.], = — y^/(z)nj3^rf-, (8.4) 

with r a contour enclosing all the ai. 

Lemma 7. Lei f be a smooth function on R satisfying If^^^'ix)] < C(l + |a;|)"'^^"' 
for 1 < n < N — 1. Then, for N > 3, there is finite constant Cn > such that 

Cm 

[ai,. . . ,aN]f < — n — rr forallai^M.. (8.5) 

l + maxi{|ai|} 

Moreover, for all 1 < n < N there is a finite constant C'j^ > such that if, 
ai < —A < for all 1 < i < n and fli > A > for all n + I < i < N , then 

C 

^ {l + X)N~2 ■ (8.6) 

Proof. We first prove (j8.6p . Using Feynman's formula 

/■ ^ 
[ai, . . . ,aAr] , = (TV - 1)! / Yjf^iCitti, . . . ,Y.f^^c^a^ 5 [l - Yj'i=iCi] W dck 
J[o,i]« L J / V ^ ^^^^^ 

= y ^ ^ f'"-'' (Etic.a.) <5 (l - Etic.) n dc, (8.7) 
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we see that [ai, . . . , un] is uniformly bounded for N > 2. Hence it suffices to 
consider the case A > 1. Similarly, a simple change of variables shows that 



[ai, . . . , unI f 

1 



N 

X 7"-^(l - jf—'S (1 - E:LiC,) S (1 - Ef=„+ic.) n ^Cfc 



fc=l 



and hence it is sufficient to consider the case ai = • • • = a„ = a and a„+i = 
aN = b. In this special case, we have 

1 f dy-W dy-' f{a)-f{b) 



where m ~ N ~ n. Note that a < —A < —1 and 6 > A > 1. The result now follows 
easily from our assumptions on /. 

Let K = maxi \ai\. To prove (|8.5I) , we may again assume that k > 1 . If max, a, — 
miui ai < k/2, then either > k/2 for all i or < —k/2 for all i, and the result 
follows from (|8.7p . If, on the other hand, max^ — min^ > k/2, the result follows 
immediately from the definition (|8.2p . with min^ in place of ai and max^ in 
place of Oat, using the boundedness of the numerator. □ 

What the lemma says is that [oi, 02, . . . , cat]/ decays at least as fast as the inverse 
of its largest argument. If all the arguments get large, but in opposite directions, 
than the decay is at least as fast as the {N — 2)th power of the inverse. 

Let now / denote the function defined in (|2.6p . It is analytic in the strip | Im z \ < 
TT. This leads to the following contour integral representation. 

Lemma 8. For R> 0, let be the contour {r — i-^, r e [— i?, i?]}U{— r + i^, r e 
[-R,R]}. Then 

fiPHA) " fWHo) = lim ^ [ f{Pz) (—^ \-] dz (8.10) 

where the limit holds in the weak sense, i.e., for expectation values with functions 
inC^{R'^)(SCf^{R'^). 

Since / is unbounded, it is important to take the difference of the two operators 
/(/?i?A) and f{f3Ho). The representation (|8.10p would not hold for the individual 
terms alone. We shall write, for simplicity, 

/(/3i?A) - fWHo) = ^ ^ fW^) (tZ^ - (8-11) 

where F is the contour z = r ± i^, r E R. This equality has to be understood as 
the weak limit 

Proof. For fixed A G M, we have 

hm -i. / /(/3z) -^dz = fm - \ , (8.12) 
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and the limit is uniform on bounded sets in A. For functions in the range of 
d{X — \Ha\) for some A > 0, this imphes that 

1 



hm 



fi(3z) (0 |(z - Ha)-'\ (^) dz = |/(/3i/A) - i| 0) 



(8.13) 



Since the operators J^^ f{l3z){z — H/^)~^dz are bounded, in absolute value, by 
C(l + |i?A|), uniformly in R, the limit extends to all functions in the form 
domain of 1 + \Ha\- The same argument applies with Hq in place of Ha- The 
form domains of 1 + I^^aI and 1 + |-ffo| are equal since A is bounded, and contain 
C;f (R'*) ® C^iR"^). This completes the proof. □ 

The diagonal entries of the 2x2 matrix- valued operator /{PHa) — /{PHq) are 
complex conjugates of each other. This follows from the unitary equivalence (|5.3p 
and the fact that /(— z) = f{z) — z. The unitary U in (|5.3p interchanges the 
diagonal entries, hence the upper left entry of /(/3-ffA) — fiPH^) equals the lower 
right entry of 



= fmA) - fimo) - P{Ha - Ho) 
= UPHa) - IWh) ~ P{Ha - Ho) 



i.U) 



Since the diagonal entries of Ha — Hq are zero, this proves the claim. In particular, 
the diagonal entries of /(/3i?A) — fiPHo) have the same trace, and hence it suffices 
to study the upper left diagonal entry, which we denote by [ • ]ii. 
Let k denote the operator 

{-ihV + hA{x))^ - ^i-h^Wix). (8.15) 

The resolvent identity and the fact that 



5:=Ha-Hq = 



A 
A 



is off-diagonal (as an operator-valued 2x2 matrix) implies that 

1 1 



Ha 



Ho 



+ h + h 



11 



where 



1 



z — k z + k 



.AT 



z — k 
1 



At- 



z - 



and 



1 



z-Ha 



1 



.At- 



1 



.At- 



1 



1 



.At- 



1 



11 



z + k z — k z + k z ~ k z + k z ~ k 
We shall first give a bound on the contribution of /3 to the integral in (|8.1ip 
Lemma 9. 



(8.16) 

(8.17) 

(8.18) 
(8.19) 

(8.20) 



f{Pz)hdz 



< Ch' 



(C)l 



..21) 
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Proof. By the triangle inequality 



f{l3z)hdz 



< / |/(/3z)|||/3||Jdz|. (8.22) 



r 



Using Holder's inequality (I3.9P and \z — Ha\ > 7r/(2/3) for z £ T, we obtain the 
bound 

II/3II1 < ^11 A||6 \\{z - krml IK^ + ^)"1L ■ (8-23) 
Moreover, since k > —ji — /i^||W^||oo, 

ioi z = r ± in / (2(3), and hence 

\\hh<j^\\A\\t. (8.25) 

It is important to get a decay faster than \z\^^ since we are integrating against the 
function | ln(l + e^^^)\ which grows linearly as z — > ~oo. We conclude that 

J^ \fWz)\ II/3II1 \dz\<C\\A\\l. (8.26) 

We can further bound ||A||6 < (27r)-'^/6||'(/;||6^^"''/'^||i||6, according to the triangle 
inequality and p.lip . Finally, HV'lle ^ C'||?/'||_f/i(c) for d < 3, by Sobolev's inequality 
for functions on the torus. This completes the proof. □ 

We continue with the bound on I2 ■ Recall the definition of the divided differences 

Lemma 10. The operator Jp f{j3z)l2dz is locally trace class, and 
-i^Tr / f{Pz)hdz- 

Pl)P2,P3e(27rZ)'' 

< Ch'-'\\t\\t\m%.^c) (8-27) 
where the function F is given by 

F{Pl,P2,P3) 

= 13 11 ^ [03,03,01,-02,-00]/ 7:7-7^ (8.28) 

and where we introduced the notation ai = l3{{q + Pi)'^ ^ '^'^'^ Pi ■— ^ 

P2 - Pa- 
One can actually show that (|8.27p holds with h^~'^ in place of h^~'^ in the error 
term. This requires substantially more effort, however, and hence we omit the proof 
for the sake of brevity. 
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Proof. The resolvent identity reads 
1 1 



{k ~ ko) 



k ' 



where 



and hence 



ka = -h^V^ 



k-ko = h'^W{x) + h'^A(x)'^ - ih^V ■ A{x) - ih'^A{x) ■ V . 



5.29) 



(8.30) 



5.31) 



We apply this to the first factor in I2 in (|8.19p . Using again the Holder inequality 
p.9|) . we can bound 



1 



1 



< 



z — k 
1 



z - ko 
[k - ko) 



1 



1 



z + k z — k z + k z — k 
\\A\\t\\iz-k)-'\\Uiz + ~kr'\\l 



.32) 



z- ko 

Using our assumptions on A and W, it is easy to see that, for z = r ± i7r/(2/3) 

for r ^ 1 



1 



- ko 
and also 



1 



ko 



-ih^VA{x) ~ ih^A{x)V) 



{h^W{x) + h^A{x)' 



<Chx 



< Ch^ X 



for r < -1 



(8.33) 



1 

T+H 



for r 3> 1 
for r <C — 1. 



As in the proof of Lemma |9l we can use (|8.24p to bound || (z — fc) ^ | 
it is not difficult to see that 



^(d-2)/(2p) 
U|-l+<i/(2p) 



for r 3> 1 
for r ^ — 1 



(8.34) 
Moreover, 

(8.35) 



for d/2 < p < 00, generalizing (|8.24p . This follows from noting that k can be 
replaced by ko for a bound, since their spectrum agrees up to o(l). With ko in 
place of fc, the result follows from evaluating the corresponding integral. In case fc 
is replaced by — fc, a similar bound holds, with the estimates for r ^ 1 and r <C — 1 
interchanged. 

For d = 3, we hence end up with a function that decays like |r|^^°/^ for neg- 
ative r and r^^ for positive r. Since we integrate against a function that decays 
exponentially for positive r and increases linearly for negative r, the remaining con- 
tour integral is finite. We conclude that ()8.32p . multiplied by f{(3z) and integrated 
over r, is bounded by C/i^~''/'^|| A||g. As in the proof of Lemma IHl we can bound 
ll^lle < C'||V'||ffi(C)'i"^^''^^ll^l|6- We have thus obtained a bound on the error made 
by replacing the first factor (z — k)~^ in I2 by (z — fco)~^. 

In exactly the same way we proceed with the remaining factors (z — k)~^ and 
(z -t- fc)~^ in /2. The only difference is that fc might now be replaced by fco in 
the factors we have already treated, but this does not affect the bounds. Also the 
number of -I- and — terms is different, but we still get a decay of at least |r|~^/'^ for 
negative r, which is sufficient for the contour integral to converge. 
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The final result is that 



1 A 1 At 1 A 1 At 1 

-A —AT —A — AT- 



z — ko z + ko z — ko z + ko z 



dz 



<oih''')Mhic) 



(8.36) 



Proceeding as above, one shows that the remaining term on the left side of (|8.36l) is 
trace class. An explicit computation of the trace per unit volume yields the desired 
result. □ 

Let us now look at the behavior of F{pi,p2,P3) for small arguments. We have 

[a,a,a,-a,~a]f = -^giia) (8.37) 
■' 16a 

with gi defined in (g^. With a ^ /3{q'^ - fi) wc thus have 



F(0,0,0) 



^ [ ^(„)4 gi(/^(g'-M)) dq 



16 



9^ - M 



(27r)° 



(8.38) 



Lemma 11. For some constant C depending on the norm oft and its derivatives 
up to order four, 



\F{pi,p2,P3) - ^^(0,0,0)1 <C{pI +pI+pI) . 



(8.39) 



Proof. We will first show that F is bounded. Using Holder's inequality, we can 
bound 

\F{p^,P2,P3)\<l3'Y[^^^^ (8.40) 

i=l 



where ^5 :— Si and 



dq 
(27r)° 



(8.41) 



for 1 < i < 4. The last factor in the integrand is bounded by a constant times 

(i + /3(<z + E;=iP0'r\ for any 1 < z < 4, as the bound (18. 5p from Lemma [7] 

shows. A further application of Holder's inequality thus shows that Si is bounded 
by Utile. 

It is thus sufficient to consider the case \pi\ < C for all i e {1, 2, 3}. We write 



tiq+p)=tiq)+p-Vt{q)+ / [p ■ V)h{q + sp){l - s)ds 



(8.42) 



and also, using partial fractions (or, equivalently, 
[l3{{q+pf ~ fi),a,b,c,d]j, 

= [f3{q^ - tJ.),a,b,c,d]j: + /3{2p- q + p^) [/3{{q + p)'^ - fi), /3{q'^ ~ fj.),a,b, c,d] . 

(8.43) 

We expand all these factors to second order in p and, in particular, iterate (j8.43p 
once more. Using the fact that the derivatives of t are in L^, as well as the decay 
bound (|8.6p of Lemma [71 we can proceed as above, using Holder's inequality, to 
arrive at the result. Note that all terms linear in p integrate to zero since t is 
reflection-symmetric. □ 
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We have 



i'{pi)ip*{p2)i'{p3)Tp*{-pi-P2-P3) ^ J \i'{x)\*dx. (8.44) 

With the aid of the Schwarz inequahty for the pi sum and Sobolev's inequahty, we 
can bound, similarly to (|6.58p - (l6.59|) . 



E 

pi,P2,P3G(2-!rZ)'' 



r{pi)r{P2H{P3m-pi-P2'P3) < c\mH-Amm ■ (8.45) 



Equal bounds hold with pi replaced by P2 and p|, respectively. Hence we arrive at 
the result 



-^Tr / f{Pz)l2dz-h''-''F{0,0,0) f \tP{x)\Ux 
2vri Jr Jc 



< Ch 



5-di\ 



(8.46) 



with F(0,0, 0) given in (|8.38p . The term involving F(0,0, 0) gives rise to the last 
term in ^J^. 

We are left with examining the contribution of Ii to the integral in (18. lip . 

Lemma 12. The operator f{l3z)Iidz is locally trace class. 

The proof of Lemma [T^ is somewhat technical and will be given in the appendix. 
Using the resolvent identity (|8.29l) we can write Ii =If + lf + Ii+lf, where 



and 



^ z - fco 



z ~ kf) z + kn z — fcn 



(8.47) 



z — ko z + ko z ~ ko 



— A^-^(ko - k)^-^A^ — 



z — kn z + ko 



2 + fco z — ko 



^ -A^^A^^^(fc-/co)- ^ 



z — ko z + ko z — ko z ^ ko 

The part consists of 6 terms involving two (fco — k) factors. Explicitly, 

z — ko z — ko z ~ ko z + ko z — ko 

+ -J—A^-iko - ~k)^^{ko ~ fc)^^At^— 

z — ko z + ko z + ko z + ko z — ko 

+ ^A^-At^-(fc - ko)^—ik ko)^— 

z — ko z + ko z — ko z — ko z — ko 

z — ko z — fco z + fco z + fco z — fco 

z — ko z — fco z + fco z — fco z — fco 

+ ^^-rr(fco - fc)-i-At^(fc - ko)^— . 

z — fco z + fco z + fc'o z — fco z — ko 



Finally, if contains all the rest. 
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A straightforward computation of the trace and the contour integral shows that 
^Tr f f{Pz)Itdz^h'-^ ^ \i^{p)\'G{hp), (8.48) 

^ pe(27rZ)'' 

where 

G[Pi) ^ P / 1 [ai,ai,~~ao\f j—^ . (8.49) 

Here, we use the notation introduced in Lemma [TUl Exphcitly, G{p) equals 
/3 f (t(g+p)+i((7))2tanh(i/3(((7+p)2-;x))+tanh(i/3(g2_^)) dq 



4 7r. 4 (q+p)2 + g2_2M (27r)'^- 

(8.50) 

Lemma 13. The function G in is bounded, twice differ entiahle at zero, and 

\G{p) - G(0) ~ i(p • V)2G(0)| < C\p\^ (8.51) 

/or some constant C depending only on J^^^ t{q)'^{l + q^)^^dq and the norm of 
t and its derivatives up to order four. 

When inserted into (|8.48l) . the term C|p|^ thus yields an error of the order 
Proof. It follows from (|8.5p of Lemma [7] that 



["^'°^'-"°]/|- l + g2+^(p + g)2 ' ^^-^2) 

hence G is bounded by Jj^at{q)'^il + q^)^^dq, which is finite by our assumption 
(|2.4p on t. To prove (|8.5ip . it thus suffices to consider the case \p\ < C. 

We expand t{q + p) up to fourth order in p. Similarly, we write ai — oq + (3(2pi ■ 
q + Pi) and expand the brackets, using the resolvent identity. The decay of the 
resulting expression for large q is governed by (j8.6p . Using in addition that t is 
reflection symmetric and satisfies (|2.3p , the result follows in a straightforward way 
from Holder's inequality. We omit the details. □ 

Recall the definition of the gi in (P?71) - (P^ . We have 

G(0) = / m^goWq' - M)) , (8.53) 

which gives rise to the expression in (j2.11l) . Using integration by parts, 
ip.VfG{0) = -^ [ tiq)[{p-V)H]iq)goif3{q'-p)) 



(27r)" 



+ ^ / {p-qrtiqfgmq'-t^)) '^^ 



+ ^ f t{qf g^q' - t^)) 7^, ■ (8.54) 

We proceed with /j" and insert ()8.3ip . We treat successively the terms propor- 
tional to W, to A^, and to A. After doing the contour integration and evaluating 
the trace per unit volume, the six terms proportional to W are 

h4-d r{piWp2)W{-pi-p2)L{hp,,hp2), (8.55) 

pi,P2 6(27i-Z)<* 



MICROSCOPIC DERIVATION OF GINZBURG-LANDAU THEORY 



37 



where 



L{puP2) = P^I Lip,,p2,q)-^ (8.56) 



and 



L{pi,p2,q) = i it{q) + t{q + pi)) {t{q + Pi) + t{q + p, + p^)) (8.57) 



Since 



we have 



X (^[ao, ao, a2, -ai] ^ + [ao, 02, 02, -ai] ^ + [aq, 02,-01, -oi] ^ j 



[a,a,a, -a]/ = ^51(0) (8.58) 



^2 2 dq 



L(0, 0) = ^ t{qY g^q' - m)) ■ (8-59) 
In the same way as in Lemma [T3l one proves that 
Lemma 14. The function L is bounded, with 

\L{puP2) - i(0, 0)1 < C [pI+pI] (8.60) 
/or some constant C > 0. 

The contribution of this term is thus 

/i^"'^L(0,0) / W{x)\ij{x)\^ dx + 0{h^-'^)\\i^\\'\j^ , (8.61) 



where we have used a simple Schwarz inequahty to bound the error term, as weU 
as the fact that W £ L^{C) and that \i}{p)\ < C||'i/;||/f2 for aU p e {2ttZ)'^. 

We obtain the same contribution with A{x)'^ in place of W{x). The terms linear 
in A are given by 

dq 

(8.62) 



1^3 d i/^*(pi)^(p2)l(-pi-p2)- / L{hpi,hp2,q)i2q + hpi + hp2) 



For small /i, we have 

/ L{hp^,hp2,q){2q + hp^ + hp2)-p-^ = hK{p2-pi) + 0{h^), (8.63) 
where IK is a symmetric d x d matrix with entries 

K,, = ii(0, 0)<5., + ^ / i(g)29,9, .g2(/3(g2 _ ^)) ^ . (8.64) 
Hence the leading term is 

h'^-'^ 2Re j ^{x)*A{x) ■ K (-^^(2;)) dx . (8.65) 
The remainder can be bounded by 

ch^-" J2 \ripimP2)Ai-pi -P2)\ (pI + pI) < ch'-''mj,^^c) (s-ee) 

pi,P2e(27rZ)'' 

using that |yl(p)||p| is summable by assumption. 
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We proceed with . First, there is the contribution of the terms linear in A in 
both factors k — ko. These we have to calculate explicitly to leading order in h, and 
then bound the remainder. They contribute 

^'"'vE / m^)\'M^)AA^)dx [ t{qfm,g,il3iQ'-f^))J^ (8-67) 

plus terms of order 0(/i^~'^)|j?/'||^i and higher. In combination, (I8.6ip for W and 
A^, (|8.67l) and (|8.54p give all the terms in (|2.12p except for the last, which came 
from F(0,0,0) in (1051) . 

For the other terms, we use (18.331) and (j8.34p . and bound all expressions involving 
z with oo-norms. We always get enough decay in either the positive or the negative 
2;-direction. Only decay in one direction is needed, as z can be replaced by — z, which 
follows from the fact that the resulting expressions are zero if fijSz) is replaced by 
f{fiz) — f{—Pz) = fiz. (Compare with the discussion at the end of the proof of 
Lemma [T2] in the appendix.) The same is true for if. We omit the details. The 
final result is then that 

^Tr ^ /(/3z) {II + If) dz = + 0{h^-%nm ■ (8-68) 

This completes the proof of Theorem [51 

9. Proof of Theorem El 

As in the proof of Theorem [21 we restrict out attention to the case d = 3, and 
leave the cases d = 1 and d = 2 to the reader. 
Since the function 

/^(^) = TTTT (9-1) 

is analytic in the strip | Im z| < tt, we can write aA with the aid of a contour integral 
representation as 

1 



OA 



2711 Jp 



z-Ha 



dz , (9.2) 



12 



where F is the contour Imz = ±7r/(2/3), and [•]i2 denotes the upper off-diagonal 
entry of an operator-valued 2x2 matrix. The integral in (|9.2I) has to be suitably 
understood as a weak limit, similarly to (|8.10p . 

Using the resolvent identity and the definitions of A and (p in (12. ip and (|2.15l) , 
respectively, we find that 

h ^ 



where 



h f ,„ . f 1 ,,,, t t ,,,,1 



^1 = 7—/ P(P^) Z — r ^o] 72 + -2 r2 [V', fco]— dz , (9.4) 



Awi Jy \z — kci ' z^ — k^ z'^ — kg ' z + k 



and 



1 



Ha 



dz , (9.6) 
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with t being short for the operator t(—ihS/). In the fohowing three lemmas we give 
bounds on these three operators. 

Lemma 15. 

\\m\\H^<ch^~'^^iA\HHc)- (9.7) 

Proof. The norm of ryi is given by 

hirm-h^-' E \^iP)\'J(hp), (9.8) 

pe(27rZ)'' 

where 

4 . od 



X |t(q)[ao, -ao, ai]p + i(p + g)[-ao, ai, -ai]p| . (9.9) 



1 2 dq 
{2tj:Y 

Here, we use again the notation [■]p for the divided differences (I8.ip - (I8.2I) . and 
ao = P{q^ — /i) and a\ = /^((g + 'pf' — /i), as in the proof of Theorem[2l Since 

. , . , p(-a)-p(a) p'(a) 1 + 2ae" - e^" 

[a, a, -a], = -[«, -a, -a], = ^ + ^ = 4«^(l + e-)^ ' ^^-^^^ 

the integrand in (j9.9p vanishes hke the fourth power of p for small p. In fact, 
using Lemma [7] and our assumptions on the regularity of t it is easy to see that 
Jij>) < C'IpI^j which yields the desired bound. □ 

We proceed with a bound on the norm of r\2- 

Lemma 16. 

\\m\\H-<Ch^-''l^mm^C)- (9.11) 
Proof. We split r]2 into three parts, 772, rf^ and 7721 where 

h'^ ( I ( 1 

'72 = -1^ / P(/3^) — r (V • ^ + ^ • V) 

27r Jp z - fco \^ z - /cq 

+ A^-(V-A + A.V) )^-dz, (9.12) 

= ^ [ pW^)^-r ( (-*v • A - 7/7^ • v)^-(fc - fco)^A 

+ (-iV • A - • V) — 5— A — ^(fco - fc) 

z — Kq z + ko 

+ A — ^-i-(-«V • A - • V) — ^(fco - k) I dz 
z + ko z + ko J z + k 

(9.13) 

and 



PiM—^ [W + m-^A - A-—{W + \A\ 



2-111 Jp z — ko \ z — k z + ko J z + k 

(9.14) 



40 



R.L. FRANK, C. HAINZL, R. SEIRINGER, AND J. P. SOLOVEJ 



We start by considering the first term in 772. The square of the norm of the 
integrand can be bounded by 



Tr 



^ ''"^\w + \A\')-l-A-^A^^iW+\An 



|z — /cop z — k |z + fc|2 2 — fc 

1 - h^W^ 



<\\A\\l\\W+\A\'"' 



\z-ko\^ 



(z-krm\\iz + krx. (9.15) 



The latter expression decays hke \z\ ^^/'^ for negative Re 2. Hence the re- 

sulting contour integral is finite, and we arrive at the bound C/i^^''/'^|j A||g < 

CH^-'^/^^/jWhiic) for the norm. 

For the second term in 773 we proceed similarly. It is important to first notice 
that 

^ 1 - (9-16) 



1 + e^^ 1 + e-l^"^ 

however, and that the 1 does not contribute anything but integrates to zero. Pro- 
ceeding as above, we conclude that 

U2\\h<Ch^-''\mlHc)- (9.17) 

The norm of 772 can be bounded in essentially the same way, and we omit the 
details. The result is that also 

\\v'2\\m<Ch'-'\mj,,^c) (9-18) 

holds 

At first sight, the norm of 772 appears to be of order h?~'^/'^ instead of h^~'^/'^. 
Note, however, that if we commute A{x) and ipi^) to the left in both terms, the 
resulting integral is zero due to the antisymmetry (I9.10p . Hence ?7f can be written 
as an integral over sums of various terms involving commutators of A(x) and ■0(x) 
with functions of — iV, and hence an additional factor of h is gained. In fact, the 
norm of 77^ can be written as 

\\V2\\m=h^-'' E ii{piW{v2)A{p^)-M{hp^,hp2,hp:i)A{-p,-p2-p3), 

Pl)P2,P3G(27rZ)'' 

(9.19) 

where M is a d x d matrix- valued function which can easily be calculated explicitly. 
An analysis as in the proof of Lemma fT5| using the antisymmetry ()9.10p . shows 
that 



n.P2.P3)\\<C{pl+pl+pl) . (9.20) 
Under our assumption that |A(p)|(l + \p\) is summable, this implies the result 

h2\\m<Ce-'"^M\mic)- (9.21) 

□ 

Finally we bound the norm of 773. 
Lemma 17. 

WmWrnl^Ch^-^Wmio- (9-22) 
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Proof. Since Ha is self-adjoint, [{z — Ha) ^]22 is bounded by 2/3/7r for z E T. Using 
Holder's inequality, we can bound the square of the norm of the integrand by 



^\\A\\l\Uz-k)-Hl\Uz + k)-Hl Vl-h^^^z-k)-' ^ . (9.23) 



This expression decays like \z\ for Rez < 0, and like \z\ ^ for Rez > 0. Hence 

1(C) 



\\V3rm<C\\A\\l<Ch^-''Ufm,c)- (9-24) 



□ 

This completes the proof of Theorem [3l 

Appendix A. Proof of Lemma [T2] 
In this appendix, we shall give the proof that the operator 

is locally trace class. As in the proof of Theorem[21 we will focus our attention on the 
case d = 3. Using the resolvent identity, we split B into three parts, B = X]i=i 
with 

and 

We shall show that each individual piece is locally trace class. 

We shall start with B2. Using Holder's inequality, the trace norm of the integrand 
is bounded by |/(/32)| times 



^ (fc - ko) 



z- ka 



WMl ll(^ - IK^ " ko)-% Wiz + (A.5) 



with 2/p + 1/g = 2/3. We shaU choose q''^ = 2/3 - e for some < e < 1/6. 
The bounds then imply that (TO]) increases as \z 

for Rez > 0, 

and decreases like Izl"^^/^"*"^ for Rez < 0. Since decays exponentially for 

Re 2 > and increases linearly for Re 2 < 0, the contour integral is absolutely 
convergent. This shows that B2 is locally trace class. In exactly the same way, one 
shows that also ^3 is locally trace class. 

For the remaining term Bi it is not possible to take the trace norm of the 
integrand and obtain a convergent integral. Hence we have to argue differently. We 
again use the resolvent identity to write 

^ ^ ^ -iko~k)^— + ^—{ko-k)^-^{ko-k)- ^ 



z + k 2 + fco 2 + fco 2 + fco z + ko z + ko z + k, 







1 (ko-k) (ko-k) (ko-k) ( A.6) 



z + ko z + ko z + ko z + k 

and, correspondingly, write Bi — B^ + B'I + Bf + Bf 
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We start with Bf. Recall that A is the sum of two terms, 

A = -^{^{x)t{-ihV) + t{-ihV)^{x)) . (A.7) 
By expanding ip in a Fourier series, we can use the triangle inequality to bound 

\\Bni<^ E mpi)mp2)\\\Dp,^pji (a.s) 

where 

Jr 2 - fco ' z + ko 

X (t(-ih\/)e''P^-'-' + e-'P=-^i(-i/iV)) -^—e'P^-'' dz . (A.9) 

z-ko 

Note that Dp-^ ^^ is translation invariant. That is, Dp^^p^ — Dp-^^p^{—ih'S/) is a 
multiplication operator in Fourier space, given by 

^Pi,P2(9) - + i(9 + hpi)) {t{q) + t{q + hp2)) 

Its local trace norm thus equals 

WDp^pJli^h-'' f \Dp,,p^{q)\j^. (A.ll) 
A simple Schwarz inequality and an application of (|8.5p show that 

independently of pi and p2. Finally, we can bound J2p l'0(p)l ^ C\\iP\\h2 for d < 3. 
We have thus shown that 




(A.13) 



This concludes the bound on the trace norm of Bf. 

The term B\ will be treated very similarly. Recall from (|8.3ip that fco — ^ is 
given by the sum of four terms. Let us first look at the two terms linear in A. We 
expand both ijj and yl in a Fourier series, and conclude that 

\\Bl\\i<^ E mPi)mP2)\\A{P3)\\\Gp,.^p...ps\\i (A.14) 

Pl,P2,P3e(27rZ)<* 

where Gp^^p2^p^{—ihS/) is a multiplication operator in Fourier space, given by 

Gpi,p2.,P3 il) ~ 

{t{q) + t{q + hpi)) {t(q - hps) + t{q + hp2 - hps)) {2q - hps) 
x[-/?(g2 _ _ hps)^ ~ ti),f3{{q + hpi)^ - M),/3((g + h{p2 ~ ps))^ ~4f 

(A.15) 
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To bound its trace norm, we can again use the bounds on the divided differences 
given in Lemma [71 We shah also need the fohowing additional bound. From the 
definition (|8.2p it follows that if / is Lipschitz continuous, then 

C / 1 1 \ ,^ , 

01,02,03,04 . < -[- r + ^ r . (A.16) 

If we apply this with oi < 02 < < 03 < 04, we obtain the upper bound (7/(6162), 
where 61 and 62 are, respectively, the largest and second largest among the numbers 
{joil, |a2|, losl, |a4|}. In combination with (j8.5|) one easily obtains the bound 



[-Piq^ - fi), -I3{{q - hp^f - m), /3((<7 + hpif - fi), (3{{q + h{p2 - Pz)f - m)] 

< I CA 17) 

- 1 + g2 1 + max{(g - /ip3)2, (g + /ipi)2, ((7 + /ip2 - /iPs)^} ■ 

A simple Schwarz inequality then shows that the trace norm of Gp-^^p^^p^ is bounded 
by C(l + |p3|) for some constant C independent of pi, p2 and P3. Since both |V'(p)| 
and |A(p)|(l + IpI) are summable, this shows that this part of is locally trace 
class. The other part, where the term quadratic in A and the term with Winko — k 
are taken into account, can be treated in the same way, using our assumption that 
|VF(p)| is summable. 

The same method can be used to show that is locally trace class. One simply 
expands all multiplication operators in Fourier series, and then bounds the resulting 
translation invariant operator. We omit the details. 

We are left with the term Bf. We first note that since 

f{Pz) = f{-/3z) + Pz (A.18) 

we can replace f{/3z) by /(— /3z) in the integrand without changing the integral, 
since the expression with /3z in place of f{f3z) integrates to zero, as can be easily 
be seen from (|8.4p . for instance. We can bound the trace norm of the integrand 
with the aid of Holder's inequality and (j8.33p - (l8.35l) by a constant times 
for Rez > 0, and \z\^/^ for Rez < 0, respectively. Since we are integrating against 
a function that decays exponentially for negative Re z and increases linearly for 
positive Rez, the integral is absolutely convergent. This proves that Bf is locally 
trace class. 
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